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INTRODUCTION 

A resurgence of intense interest, research and applica- 
tions mark the beginning of what may be considered the second 
generation of filamentary composite materials. Such interest 
and activity are as well-founded now as they were at the out- 
set of the composites era more than twenty years ago. The 
possibility of using relatively brittle materials with high 
modulus, high strength, but low density in composites with 
good durability and high tolerance to damage and which, when 
they do fail, do so in a noncatastrophic manner, has been 
shown feasible, and the full potential is only just begin- 
ning to be realized. 

The promise of substantially improved performance and 
potentially lower costs provides the driving force behind 
continued research into fiber reinforced composite materials 
for application in aerospace hardware. 

Much progress has been achieved since the initial devel- 
opments in the mid 1960's. Applications to primary structure 
have been rather limited on operational vehicles, mainly be- 
ing utilized in a material-substitution mode on military air- 
craft. More extensive experiments, as a part of NASA's in- 
fluential ACEE program, are currently underway on large air- 
planes in commercial passenger operation and in a few military 
developments, such as the AV-8B which has seen only limited 
service use and the X-29 which is in its initial flight 
testing program. 
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A strong technology base is required to fully exploit 
composites in sophisticated aerospace structures. NASA and 
AFOSR have supported expansion and strengthening of the tech- 
nology base through programs which advance fundamental know- 
ledge and the means by which this knowledge can be success- 
fully applied in design and manufacture. 

As the technology of composite materials and structures 
moves toward fuller adoption into aerospace structures, some 
of the problems of an earlier era are being solved, others 
which seemed important are being put into perspective as 
relatively minor, and still others unanticipated or put aside 
are emerging as of high priority. The purpose of the RPI 
program as funded by NASA and AFOSR has been to develop cri- 
tical advanced technology in the areas of physical properties, 
structural concepts and analysis, manufacturing, reliability 
and life prediction. 

Our approach to accomplishing these goals is through an 
interdisciplinary program, unusual in at least two important 
aspects for a university. First, the nature of the research 
is comprehensive - from fiber and matrix constituent proper- 
ties research, through the integration of constiuents into 
composite materials and their characterization, the behavior 
of composites as they are used in generic structural compo- 
nents, their nondestructuve and proof testing and, where the 
state of the art will be advanced by doing so, extending the 
research effort into simulated service use so that the com- 
posite structure's long-term integrity under conditions 
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pertinent to such use can be assessed. Inherent in the RPI 
program is the motivation which basic research into the struc- 
tural aspects provides for research at the materials level, 
and vice versa. 

Second, interactions among faculty contributing to pro- 
gram objectives is on a day to day basis without regard to 
organizational lines. These contributors are a group wider 
than that supported under the project. Program management 
is largely at the working level, and administrative, scien- 
tific and technical decisions are made, for the most part, 
independent of considerations normally associated with aca- 
demic departments. Depending on the flow of the research, 
this kind of involvement includes faculty, staff and students 
from chemistry, civil engineering, materials engineering, 
aeronautical engineering, mechanical engineering and mechanics. 

A third aspect is a newly instituted initiative to in- 
crease involvement of NASA's Research Center scientists and 
engineers in the program at RPI, and vice versa. This has 
required, first, identification of individual researchers 
within NASA centers whose areas of interest, specialization 
and active investigation are in some way related to those of 
RPI faculty supported under the subject grant. Second, a 
program of active interchange has been encouraged and the 
means by which such interaction can be fostered is being 
sought. Important benefits envisioned from this increased 
communication include a clearer window to directions in aca- 
demia for NASA researchers; opportunities to profit from NASA 
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experience, expertise and facilities for the faculty so in- 
volved and an additional channel for cross-fertilization 
across NASA Research Center missions through the campus pro- 
gram. 

Overall program emphasis is on basic, long-term research 
in the following categories: (a) constituent materials, (b) 

composite materials, (c) generic structural elements, (d) 
processing science technology and (e) maintaining long-term 
structural integrity. Depending on the status of composite 
materials and structures research objectives, emphasis can 
be expected to shift from one time period to another among 
these areas . Progress in the program will be reported in the 
following pages under these headings. Those computer method- 
ology developments are also undertaken which both support RPI 
projects in composite materials and structures research in 
the areas listed above and which also represent research with 
the potential of widely useful results in their own right. 

In short, the NASA/AFOSR Composites Aircraft Program is 
a multi-faceted program planned and managed so that scientists 
and engineers in a number of pertinent disciplines at RPI will 
interact, both among themselves and with counterpart NASA 
Center researchers, to achieve its goals. Research in basic 
composition, characteristics and processing science of compos- 
ite materials and their constituents is balanced against the 
mechanics, conceptual design, fabrication and testing of 
generic structural elements typical of aerospace vehicles so 
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as to encourage the discovery of unusual solutions to present 
and future problems. In the following sections, more detailed 
descriptions of the progress achieved in the various component 
parts of this comprehensive program are presented. 
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PART II 


II— A MECHANICAL PROPERTIES OF HIGH PERFORMANCE CARBON FIBERS 
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II-A MECHANICAL PROPERTIES OF HIGH PERFORMANCE CARBON FIBERS 

Senior Investigator: R. J. Diefendorf 

1. Introduction 

The studies of modulus distribution, strength distribu- 
tion and residual stress measurement within a single fiber 
are of great interest. They may lead to the explanation of 
the anomalous behavior of modulus and strength with increasing 
heat-treatment temperature. Also, the results may hold the 
key to the realization of simultaneously developing high 
modulus and high strength fibers. 

2. Status 

Theoretical equations for predicting residual stresses 
in radial, hoop and longitudinal directions were derived from 
linear elasticity based on the perfect onion skin model and 

m * 

recorded in a recent semi-annual progress report . Resid- 
ual stress measurements in high modulus, PAN-based carbon 

fiber, as well as modulus distribution in Type I, Type II and 

n 2 3i 

Type A PAN-based carbon fibers, were explored and reported ' ’ 
Strength distributions within a single fiber of Type I, Type 
II and Type A PAN-based carbon fibers were established. Flaws 
on the fracture surface of a single fiber also were identi- 
fied. 


Numbers in brackets in this section refer to the references 
which are listed on page 23. 
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3. Progress During Report Period 

Hercules HMS, HTS , AS-4 and AS fibers and Union Carbide 

T-300 fiber were used in this study. The surface of the 

[2 3] 

carbon fiber was removed by electrochemical etching ' . 

The diameters of the fiber, both before and after etching, 

T41 

were measured by a laser diffraction technique . The 
strength of single filaments was measured at a gauge length 
of 25.4mm under test conditions that allowed recovery of the 
primary fracture surfaces needed for identifying critical 
flaws . 

fl 3i 

The surface layer of carbon fiber has higher modulus ' 
and, hence, will carry a higher fraction of load than the 
rest of the layers. Therefore, by removing the outside layer, 
the strength will also change. This strength distribution 
will vary among types of fiber and depends on the magnitude 
of local modulus, residual stress and flaw distribution. 

The strength distributions of HMS, HTS, AS-4, AS and T- 
300 fibers are shown in Figures II-A-la through II-A-le and are 
based on the average value in each radius ratio group (via, 
group 1 for as-received fibers r/r Q = 1, group 2 for etched 
fibers with radius ratios ranging from 0.95 to 0.99, then 
r/r Q = 0.9-0.94 for group 3 and so on). By performing statis- 
tical analyses, such as Bartlett's test, ANOVA, Dunnett's 
test and the t-test^, the group mean values were properly 
compared with each other. For HMS fibers, the mean strength 
of group 1 (as-received) is higher than those after etching. 




NORMALIZED RADIUS <r/r 0 > 


Figure II-A-la. Strength Distribution as a Function of Normalized Fiber Radius for HMS Fibers. 
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Figure II-A-lb. Strength Distribution as a Function of Normalized Fiber Radius for HTS Fibers. 




Figure II-A-lc. Strength Distribution as a Function of Normalized Fiber Radius for AS-4 Fibers. 
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Figure II-A-ld. Strength Distribution as a Function of Normalized Fiber Radius for T-300 Fibers. 
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Figure II-A-le. Strength Distribution as a Function of Normalized Fiber Radius for AS Fibers. 
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This finding is related to the structure of HMS fiber and the 
effect of residual stress distribution within the fiber. Be- 
cause the layers of higher local modulus were etched away 
gradually, the average modulus of the remaining fiber became 
smaller and smaller. Simultaneously, the residual stress is 
relieved. Therefore, the surface compressive stress de- 
creased. In terms of a flaw mechanism, most of the fracture 
surfaces of this high modulus species showed gross interior 
voids (see Figure II-A-2) , and fracture is believed to orig- 
inate from these voids. By removing the outer layer, the 
flaw originally in the interior was exposed as "surface" flaw 
(see Figure II-A-3) . Also, less compressive stress protects 
this "surface" flaw from propagating, therefore, the strength 
drops. By comparing the mean strengths, the group pairs 
(2,3), (3,4), (4,5) and (5,6) showed that they share the same 
mean value. That is, the strength of HMS fiber dropped after 
the surface was removed and then stayed constant (as shown 
in Figure II-A-la) . 

HTS carbon fibers, as shown before^, have less of a 
modulus gradient and lower HTT than that of HMS fibers; there- 
fore, less residual stress is expected. The strength distri- 
butions of HTS fibers are shown in Figure II-A-lb. Although 
the nominal mean values of strength in etched fiber groups 
are higher than those of as-received fibers, they are not sig- 
nificantly different statistically (see Figure II-A-lb) . The 
fracture surfaces of HTS as-received fiber (Figure II-A-4) 
showed that failure initiated from a surface flaw. No serious 





Figure II-A-3. SEM Fractography of HMS 
(Etched) Carbon Fiber. [shi 
that the flaw originally in 
the interior was exposed as 
"surface" flaw) 
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Figure II-A-4. The Failure Initiated from Surface Flaw 
for HTS (As-Received) Carbon Fibers. 
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interior voids, as in HMS fibers, are observed. This implies 
that less residual stress is present, otherwise, the surface 
compressive residual stress tends to inhibit the propagation 
of the crack from the surface. After the surface layers are 
removed, the average modulus is slightly decreased. There- 
fore, for constant strain the strength decreases only slight- 
ly. Since surface flaws may be removed when etching, the 
strength might be higher than, or stay at the same value, as 
group 1 fibers. 

The past report ^ indicated that the modulus gradient 
of Type A (AS-4, AS and T-300) fibers is among the lowest of 
the carbon fiber types. Also, as the heat treatment temper- 
ature of Type A fibers is only 1100 °C, the effect of residual 
stress is expected to be low. The group mean strength of 
AS-4 is all equal within the groups (as shown in Figure II- 
A-lc) . The grand mean strength was 504 Ksi for a radius 
ratio down to .65 and with a pooled standard deviation of 
150 Ksi. 

The same results as with AS-4 fibers are found in T-300 
fibers by using the Barlett's test and ANOVA (see Figure II- 
A-ld) . However, the strength of AS fibers did not follow the 
same trend as other Type A fibers such as AS-4 and T-300 (see 
Figure II-A-le) . The mean strength of AS as-received fibers 
is only 355 Ksi, which is quite low compared to manufacturer's 
data (410 Ksi). Hence, the distribution could be due to scat- 
tering of data. 
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In general, there are two distinct types of distribu- 
tions of strength among types of carbon fibers. Type I fi- 
bers have a strength gradient when etched, while Type II and 
Type A strength gradients were not observed. These strength 
distributions are related to the microstructure of the fi- 
bers . For Type I with a skin/core structure , the strength 
is higher at the surface; then it degrades toward the core, 
as the surface layer is removed. In addition, due to the 
relief of residual stress and the exposure of the interior 
flaws to the surface after etching, the strength of the 
etched fiber further decreases . Interior voids are found in 
Type I which are the main failure initiator for Type I. For 
Type II and Type A, no significant variations were observed 
in fiber strength as the surface was removed. Surface flaws 
are the main failure initiator for Type II and Type A fibers. 
Finally, no serious interior voids were found either in Type 
II or Type A fibers. 

4. Plans for Upcoming Period 

During the next reporting period, effort will concen- 
trate on the following: 

a) finding a suitable technique for uniformly etch- 
ing pitch based carbon fiber as well as for re- 
sidual stress measurement and 

b) performing a fundamental study of the electro- 
chemical etching mechanisms of carbon fiber by 
using a highly annealed pyrolytic graphite (HAPG) 
as the sample and methanol as the electrolyte. 
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III— A FATIGUE IN COMPOSITE MATERIALS 

Senior Investigator: E. Krempl 

1. Introduction 

The deformation and failure behavior of graphite/epoxy 
tubes under biaxial (axial and torsion) loading is being in- 
vestigated. The aim of this research is to increase basic 
understanding of and provide design information for the bi- 
axial response of graphite/epoxy composites. 

2. Status 

[ 1 ] * 

In the previous progress report phenomenological 

damage laws based on residual strength or residual life were 
discussed, and a model was proposed from which various pre- 
viously published models could be recovered by specialization. 
Specimen fabrication and tensile as well as fatigue testing 

of Gr/E [±45] tubes, in support of the theory, comprised 
s 

most of the activity in the reporting period. 

3. Progress During Report Period 

The specimen design and the [±45] lay-up developed at 

s 

[21 

Rensselaer Polytechnic Institute were used to fabricate 
about fifty thin-walled tubes of graphite/epoxy. The details 
of this procedure are given in Reference [2] or in previous 
progress reports. 


Numbers in' brackets in this section refer to the references 
which are listed on page 31. 
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Subsequently, static tests and fatigue tests were per- 
formed following the procedures of Reference [2] except that 
most of the tests were performed without extensometry . 

The damage evolution equation proposed in Reference [1] 
is 

^=g(D)f(s) (1) 

where: 

D is the dimensionless damage, 

N is the number of cycles, 

s is the dimensionless stress amplitude and 
g and f are positive functions. 

Equation (1) has no repository for mean stress and there- 
fore cannot reproduce such effects. Accordingly, tests were 
run to check on the influence of mean stress. Figure III-A-1 
shows that there is a small influence of mean stress for R 

values (R = a . /a where a and a . are the maximum and 
mm max max mm 

minimum stresses of a cycle, respectively) ranging from R = -1 
(completely reversed), R = 0- (fluctuating compression) to R = 0 
(fluctuating tension). The normalizing factor of 145 MPa is 
the average tensile strength obtained with six specimens. 

The following tests were all performed at R = 0 and a 
frequency of 5 Hz. First, the drop of residual strength with 
cycles was determined with the results given in Figure III- 
A-2. The circles indicate the residual strength after a fa- 
tigue test with s = 32% which was terminated at the indicated 
number of cycles. Subsequently, a tensile test was performed, 
and its result is plotted on the ordinate. 



NORMALIZED STRESS AMPLITUDE 


in 



Figure 


III-A-1 


The Influence of Mean. Stress on the Axial Fatigue Behavior of [±45]^ Tubes. 

( R = CT min/ a max where a and a • are the maximum and. minimum stress of a cycle, 
respectively^ 34 s = stress amplitude/145 MPa) 
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Figure III-A-2. Residual Strength (made dimensionless by dividing 
by 145 MPa) After Prior Uniaxial Fatigue Cycling 
at R = 0 with s = 32%. (Dashed line represents 
the maximum stress in a fatigue cycle.) 
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Figure III-A-2 shows the drop off of residual strength 
with cycles. The dashed line is the maximum fatigue stress 
of the fatigue test. At the cycles-to-f ailure corresponding 
to this maximum stress level, the residual strength is ex- 
pected to be equal to it. The average life at an amplitude 
of 32% is 11,000 cycles. From this information and from the 
results shown in Figure III-A-2, it must be concluded that 
the residual strength decreases rapidly at the end of the 
life. 

The results of Figure III-A-3 are obtained with speci- 
mens subjected to 5000 cycles at R = 0 with s = 32% and sub- 
sequent combined loading. The triangles represent the com- 
bined residual strength of the damaged specimens. The circles 
are results taken from [1] and represent the combined strength 
of virgin specimens. 

It is tentatively concluded that the strength reduction 
due to uniaxial tension-tension fatigue cycling is uniform 
and affects the fracture strength in the first quadrant of 
the axial and shear stress space equally. 

4. Plans for Upcoming Period 

The experimental investigation outlined above will be 
continued using other test conditions. 

5. References 

1. 46th Semi-Annual Progress Report, Composite Structural 
Materials, Rensselaer Polytechnic Institute, Troy, NY, 
August 1984, pp. 29-33. 
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Figure III-A-3. Biaxial Fracture Strength of Virgin, 0, and Pre 
damaged. A, Specimens. (The damage was introduced 
during 5000 prior cycles at R = 0 and s = 32%. Re 
suits of virgin specimens are from [l] . a = axial 
stress, t = shear stress.) 
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2. T. M. Niu, "Biaxial Fatigue of Graphite/Epoxy [±45] g 

Tubes", Dr. Eng. Thesis, .Rensselaer Polytechnic Institute 
Troy, NY, May 1983. 


6 . Current Publications or Presentations by 
Professor Krempl on this Subject 

"Biaxial Fatigue Properties of Thin-Walled Composite Tubes", 
with D. M. Elzey, T. Ayar and R. G. Loewy. 

To be presented at the Helicopter Fatigue Special- 
ists' Meeting, American Helicopter Society Midwest 
Region, St. Louis, MO, October 16-18, 1984. 

"Time-Dependent Deformation and Fatigue Behavior of [±45] g 
Graphite-Epoxy Tubes under Combined Loading", with T.-M. Niu 

To be presented at the ASTM Symposium on Composite 
Materials: Fatigue and Fracture, Dallas/Ft. Worth, 

TX, October 24-25, 1984. 

"Biaxial Fatigue and Damage Accumulation in Thin-Walled Gra- 
phite/Epoxy Tubes", with D. An. 

To be presented at the IUTAM Symposium on the Me- 
chanics of Damage and Fatigue, Technion-Israel In- 
stitute of Technology and Tel Aviv University, 

July 1-4, 1985. 
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III-B EXPERIMENTAL AND THEORETICAL STUDIES OF MOISTURE AND 
TEMPERATURE EFFECTS ON THE MECHANICAL PROPERTIES OF 
GRAPHITE/EPOXY LAMINATES AND NEAT RESINS 

Senior Investigator: S. S. Sternstein 

1. Introduction 

During the present report period a new series of mechan- 
ical tests were initiated, as described herein. 

Previous reports have described our small strain visco- 
elastic studies on carbon/epoxy laminates. Techniques such 
as stress relaxation, creep and dynamic mechanical spectros- 
copy were used to investigate the in situ behavior of the 
epoxy as a function of temperature and moisture history. One 
of our findings was that the so-called loss factor (i.e., tan 
6) or ratio of out-of-phase modulus to in-phase modulus was 
sensitive to damage. For linear viscoelastic behavior it can 
be shown that loss factor is equal to the ratio of energy dis- 
sipation to energy storage per quarter cycle of sinewave. 

For nonlinear viscoelastic behavior, or for dissipative 
processes which do not have fading memory or obey time-trans- 
lation superposition, it is difficult or impossible to assign 
any specific meaning to the loss factor, as defined by phase 
angle measurements. How then does one characterize high 
load, large deformation, nonrecoverable dissipative phenomena? 

One technique which directly measures energy dissipation 
is the cyclic hysteresis test. In this test, the sample is 
repeatedly loaded and unloaded between specified limits with 
continuous recording of load versus displacement. The area 
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inside the loop so formed is hysteresis or energy dissipated. 
The magnitude of the hysteresis is expected to vary with the 
number of cycles performed, the cycle limits and the previous 
history of the sample. 

If the process being observed is linear viscoelastic, 
then the energy dissipation will vary as the square of the 
limit stress (or strain) . On the other hand, phenomena such 
as crack formation in the matrix, fiber failure, fiber reor- 
ientation, interfacial failure, etc. , are not viscoelastic 
(let alone linear viscoelastic). Thus, it is not clear as 
to how hysteresis should vary with the magnitude of stress 
or strain. However, one can anticipate that increasing dam- 
age will affect hysteresis and the rate of change of hyster- 
esis with stress (or strain level) . 

2. Status 

A preliminary set of experiments have been performed 

/ 

which appear to be quite promising. We have used the cyclic 
hysteresis method to measure energy dissipation on two dis- 
tinctly different samples, one being a thermoset matrix com- 
posite (carbon-epoxy) , the other a thermoplastic matrix 
(polysulfone-carbon) composite. Both sets of samples were 
prepared with a [±45]^ s stacking sequence and tested in ten- 
sion. The approximate sample dimensions were: gage length = 

1.4", width = 0 . 75 " , thickness = 0 . 065" . The samples were 
tested using a constant crosshead speed of 0.2 cm/min. which 
gave a nominal strain rate of 5%/min. Compound limits were 
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used for this cycle, that is, the lower cycle point was set 
for zero load while the upper cycle point was set for a con- 
stant displacement. Cycling was performed for each sample 
until a stable (repetitive) hysteresis loop was obtained. A 
new sample specimen was then used for the next level of up- 
per displacement limit. 

3. Progress During Report Period 

Typical results are shown in Figure III-B-1 for the 
thermoset matrix composite and in Figure III-B-2 for the 
thermoplastic matrix composite. The first cycle hysteresis 
loop is shown as well as the final repetitive (stable) hys- 
teresis loop. After reaching the stable hysteresis loop 
condition, the samples were annealed in an oven at 100°C 
for twenty-four hours and retested. These data are shown 
in Figures III-B-1 and 2 as "recovery". It is interesting 
to note that the thermoplastic matrix composite (PS) shows 
a stiffening of the recovery curve relative to the original 
curve whereas the thermoset composite (Fiberite) shows a 
permanent softening. 

The areas inside the hysteresis loops were measured and 
normalized by the volume of the sample, thereby obtaining 
energy dissipation per unit volume (in units of lb-in/in 3 or 
psi) . The results are shown in Figure III-B-3 where energy 
dissipation is plotted versus the peak stress at the upper 
cycle point. It is interesting to note that the thermoset 
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Figure III-B-1. Typical Hysteresis Loops Showing Load Versus Displace 
ment for Fiberite, a Carbon-Epoxy (Thermoset) Laminate. 
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Figure III-B-2. Typical Hysteresis Loops Showing Load Versus Displace 
ment for a Carbon-Polysulfone (Thermoplastic) Laminate. 
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Figure III-B-3. Normalized Energy Dissipation Versus Peak Stress 
Level in the Cyclic Hysteresis Test. 
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(Fiberite) shows more nonlinearity with stress level than 
does the thermoplastic. It remains to be determined whether 
or not this represents an increased rate of damage produc- 
tion. Figure III-B-3 shows both the first cycle and stable 
cycle energy dissipation values. 

The slopes of the curves in Figure III-B-3 are equal to 
or greater than four. Thus, energy dissipation is varying 
with the fourth power of stress. This is in marked contrast 
with the slope of two expected for linear behavior. We are 
currently continuing these tests at lower stress levels to 
determine if the curves of Figure III-B-3 bend over and ap- 
proach a slope of two at lower stress values. If this hap- 
pens, it may be possible to isolate the stress level at which 
damage processes initiate. 

The energy dissipation is also shown versus upper cycle 
displacement (strain) in Figure III-B-4. Curiously, the 
thermoset and thermoplastic composites show surprisingly sim- 
ilar behavior when compared on a peak strain basis. 

4. Plans for Upcoming Period 
During the next period we plan to continue this work 
and extend it to several other thermoplastic matrix composites 

5. Current Publications or Presentations by 
Professor Sternstein on this Subject 

"Mechanical Characterization of Composites" 

Presented the invited lecture at Composites and 
Paper Technology Interchange, Stockholm, Sweden, 

June 25-29, 1984. 
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Figure III-B-4. Normalized Energy Dissipation Versus Peak Strain 
Level in the Cyclic Hysteresis Test. 



"Mechanical Characterization and Failure of Composites" 

Presented as an invited lecture at the Materials 
Research Council Meeting, La Jolla, CA, July 19-20, 
1984. 

"Viscoelastic Characterization of Composites" 

Presented at the 188th National ACS Meeting, Phila- 
delphia, PA, August 26-31, 1984. 



44 



45 


III-C NUMERICAL INVESTIGATION OF THE MICROMECHANICS OF COM- 
POSITE FRACTURE 

Senior Investigator: M. S. Shephard 

1. Introduction 

To understand mechanisms of failure in composites it is 
necessary to develop insight into their micromechanical be- 
havior, including interactions between matrix and fibers as 
the load is increased from zero to that corresponding to 
failure. Investigation of these phenomena, either experi- 
mental or numerical, is difficult. The purpose of this 
project, carried out primarily by Nabil Yehia, is to develop 
a finite element program capable of performing the numerical 
investigations of two-dimensional micromechanical failure of 
composites. 


2. Status 

A generalized program for the two-dimensional analysis 
of static crack growth problems has been developed and tested. 
In addition to an examination and development of fracture cri- 
teria needed to track cracks in composites at the micromechan- 
ical level, the program employs fully automatic finite element 
mesh generation. This is an entirely n ew approach in program 
structure that was taken to allow for the automatic tracking 
of discrete crack growth. 
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3. Progress During Report Period 
During the last reporting period, our effort concen- 
trated on (a) completing the integration of the fully auto- 
matic mesh algorithms and (b) running of test cases to check 
the program's capabilities and to compare results. 


a. Integration of Automatic Crack Tracking System 

Previous progress reports have discussed the functions 

of the automatic crack tracking system which include: 

[ 11 * 

1) a geometric based program control structure , 


2) automatic mesh generation accounting for the 

[ 1 ] 


[i-3i 

presence of crack tips and 


3) finite element based fracture criteria 
The geometric based program control structure represents the 
link to the various parts of the program. It tells the mesh 
generator how and what to mesh and how to relate the loads, 
material properties and boundary conditions to that mesh. 

The analysis and propagation results are tied back to geo- 
metry which is updated using the geometry modification oper- 
ators outlined in the last progress report. Because these 
operators must account for general geometric modifications, 
an efficient and flexible data structure must be employed. 

Within the program geometry, entities are stored in 
linked list structures^ 1 ] This data structure can account 
for geometric modification by simply adding data and shifting 
pointers within the list. For example. Figure III-C-1 


Numbers in brackets in this section refer to the references 
which are listed on page 65. 
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Point Structure 


Curve Structure 




Figure III-C-1. Data Manipulation for Dividing a Curve into Two Curves. 
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demonstrates the changes to the point and curve data struc- 
ture that occur when a curve is divided into curves. In this 
case, an entry is added to the curve list, and pointers in 
the point list are altered. No data needs to be shifted in 
this process. 

In addition to the geometric operators needed to track 
crack propagation, features were needed to account for modi- 
fications to attribute (load, material property and boundary 
condition) information and to insure that a mesh would be 
generated properly for the given geometry^ 1 ] Specific fea- 
tures included an ability to: 

1) automatically change mesh size as a crack tip 
approaches a boundary, 

2) eliminate any crack tip that hits an exteriod 
boundary and 

3) split distributed loads and boundary conditions 
on curves cut by crack tips. 

The final step in integrating the program was incorpor- 
ating the crack nucleation and propagation algorithm that 
ties the finite element analysis and fracture prediction back 
to the geometry. 

Once the user has defined the geometry of his problem 
and all the associated attributes related to geometry, the 
automated analysis is entered and proceeds as follows : 

Step 1 : 

• Read the problem definition and all attributes from the 

binary files to the COMMON blocks. 
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Step 2: 

• Generate the finite element mesh and minimize the ele- 
ment wave front. 

Step 3: 

• Formulate the element stiffness matrices and solve them. 
Step 4: 

• Evaluate element stresses at the numerical integration 
points. 

Step 5: 

• If requested, enter the following crack nucleation 
routine: 

If maximum principal stress is less than the ten- 
sile strength, continue, else compute nucleation 
scaling factor F = o t /a 1 . If F is greater than 
the previous one, Fp , continue, else set F = F p 
then update the nucleation point coordinates and 
continue. 

• Else continue. 

Step 6 : 

• Set the nucleation scaling factor F nu = F and print out 
the nodal displacements and element stresses. 

Step 7 : 

• If crack tips exist, enter the following crack propaga- 
tion routine: 

a. Calculate the stress intensity factors at each 
crack tip. 

b. Calculate the crack propagation angle 0 . 

c. Calculate the crack propagation increment 
r b (0 o ) in the modified T-criterion and/or r in 
the S-criterion, in the predicted direction of 
propagation 0 O . Then, compare r b or r with 
the material critical value, r cr , and evaluate 
the propagation scaling factor, F, at each 
crack tip as; 

-F = F/r cr /r b for the modified T-criterion or 

-F = /? 7r for the S-criterion. 
cr 

Set the propatation scaling factor for the cur- 
rent iteration, F pr , equal to the smallest F. 

• Else continue. 
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Step 10: 

Continue with Step 2 . 
b. Results of Examples ■ 

The previous progress report presented two example prob- 
lems. Four additional example problems^ 1 -* are presented 
here. 

Example 1: 

In this example, the slant crack problem is con- 
sidered under remote, uniform, uniaxial tension. 

The geometry considered has the dimensions 

0.3*0.09*0.0015 in. which are the same as those of 

[4] 

the specimens tested by Theocaris et al. A 

crack with total length 2 a = 0.016 in., making an 
angle of 40 degrees with the load axis, is con- 
sidered. Since no material properties for the test 
specimens have been reported by Theocaris et al., 
the following properties are assumed, E = 250 psi, 
a^ = lpsi, K lc = 1 psi/in. and v = 0.3. The geome- 
try and the finite element idealization of this 
problem are shown in Figure III-C-2. The initial 
crack propagation angle predicted by the modified 
T-criterion with the quarter point displacement 

technique is found to be -62.77° compared to 59° 

r 4 1 

reported experimentally by Theocaris et al . 

(i.e., 6% error). The S-criterion yields an angle 
of -51.48° (i.e., -12.7% error). With the modified 
T-criterion, the fracture load normalized with mode 
I fracture load is found to be 1.443 compared to 
1.098 reported experimentally in Reference [4] 

(i.e., 31% error). The S-criterion yields a nor- 
malized fracture load of 1.548 (i.e., 41% error). 

The propagation was always unstable. Figure III-C-3 






Figure III-C-3. Crack Propagation Path as Predicted by the Pro- 
gram CPFFEP in Comparison with the Path Experiment 
ally Determined by Theocaris et al. 


shows the propagation path predicted by both the 

modified T-criterion and the S-criterion and the 

approximate propagation path determined experi- 
[4] 

mentally . 

Example 2 : 

In this example, the slant crack geometry experi- 

[51 

mentally tested by Sih is considered under plane 
stress conditions. The geometry is 18*9*0.1875 in. 
with a total crack length of 2a = 2 in. making a 40 
degree angle with the load axis. The material 
properties are E = 100ksi, a y =29.2ksi, K lc = 29.2 
ksi/in. and v = 0.333. The geometry and finite el- 
ement idealization of this problem are shown in 
Figure III-C-4. The initial crack propagation local 
angle predicted by the modified T-criterion using 
the quarter point displacement technique is found to 
be -62.83° compared to -55.6° reported experimental- 

[51 

ly by Sih (i.e., +13% error). The S-critenon 
yields a normalized fracture load of 1.48 (i.e., 
+15.6% error). Note that better agreement between 
the numerical predictions and the experimental re- 
sults is observed in this example. The propagation 
was always unstable. Figure III-C-5 shows the pro- 
pagation path of the upper crack tip predicted by 
both the criteria and initial propagation path de- 

[51 

termxned experimentally ; 

Example 3 : 

This example is the 1 in. thick plate with two holes 
as shown in Figure III-C-6. The starting mesh for 
this example is shown in Figure III-C-7. As indi- 
cated in Figure III-C-6, the plate is put into ten- 
sion by a uniform load of 1 lb/in. , and the circular 
hole has a small crack emanating from it. A plane 
stress condition is assumed with the following 
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Figure III-C-4. Slant Crack Geometry Tested by Sih 

and Its Finite Element Mesh. 
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Figure III-C-7. 


The Starting Mesh for the Plate with Two Holes. 
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material properties, E = a = 1 psi, K =1 psi/in. 

JL J- C 

and v = 0.333. The meshes for two steps of crack 
propagation are shown in Figure III-C-8. The frac- 
ture process, with a total of six propagation in- 
crements, as shown in Figure III-C-9, has been 
tracked using both the modified T-criterion (MT) 
and the S-criterion (S) with stress intensity fac- 
tors calculated using the quarter point displace- 
ment technique. The predicted fracture loads for 
the next propagation increment are found to be p = 
7.23, 9.66, 4.59,1.27, 0.66 and 0.49 lb/in. Thus , 
the crack growth is stable for the first two incre- 
ments, as indicated by positive load increments 
predicted for the next propagation. The rest of 
the propagation increments are unstable as indicated 
by the negative load increments. For this example, 
the crack propagation load and path predicted by the 
modified T-criterion are very close to the ones pre- 
dicted by the S-criterion. However, when employing 
the assumption of [5] for choosing the positive 
angle of 0 O for cracks under compression (i.e. , 
negative K I and K ) , the first increment predicted 
by the S-criterion is found to occur at a local tip 
angle of 179° (i.e. , the crack reverses its direc- 
tion) . This is 180° different from the propagation 
angle predicted by the modified T-criterion. This 
difference vanishes if the restriction of [5] is 
removed, which yields the propagation path predicted 
by both criteria (see Figure III-C-9). 

Example 4 : 

In this example, a single fiber-matrix composite 
model is considered. It consists of a 2 in . diameter 
fiber embedded in an annular sheath of a 4 in. thick 
matrix with a 2 in. interface crack and loaded in the 



a) Second Propagation Increment 



Figure III-C-8. Finite Element Mesh of the Plate. 
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Figure III-C-9. Crack Propagation History Predicted by the Modified 
T-Criterion and the S-Criterion of the Plate 
with Two Holes. 
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fiber's axial direction. To simulate this axisym- 
metric problem, a two-dimensional plane strain anal- 
ysis of one quarter of the model is considered. 
Figure III-C-10 shows the model geometry and boun- 
dary conditions. The initial finite element mesh 
is shown in Figure III-C-11. For fracture analysis 
of this problem the modified T-criterion for inter- 
face crack propagation^ is employed, with the 
stress intensity factors being calculated by the 
crack flank displacement technique^] 

Two composite materials are considered, E f /E m = 2 
(E is the Young's modulus) which is almost the case 
of a ceramic/ceramic composite, and E f /E m = 100 which 
is similar to the case of a high modulus graphite/ 
epoxy composite. The Young's modulus of the matrix 
E m is assumed to be 1 psi while the Poisson's ratio 
v m = v f = 0 . 3 is considered in the two cases. 

Two loading conditions are considered, constant 
longitudinal strain (i.e., e ma trix = e fiber^ and un i~ 
form uniaxial load, both applied remote from the 
crack (see Figure III-C-12) . Under the constant 
longitudinal remote strain condition, crack closure 
occurred, leading to a very high load scaling fac- 
tor required for crack propagation along the inter- 
face. In the case of uniform uniaxial load, four 
subsequent equal propagation increments are pre- 
dicted at p = 15, 11.58, 9,13 and 7.34 lb/in for the 
case of E f /E m = 2 and p = 2.38, 1.73, 1.32 and 1.06 
lb/in for the case of E f /E m = 100. Therefore, in 
this model, the resistance to interface crack propa- 
gation increases as the Young's modulus of the ma- 
trix approaches that of the fiber. Although the 
propagation load varies substantially between the 
two load cases in the two composites considered, 




Figure III-C-10. A Single Fiber-Matrix Composite Model 
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Figure III-C-11. 


The Initial Mesh for the Single Fiber-Matrix 
Composite. 





65 


the propagation has been always found unstable. 

The meshes for the second and fourth propagation 
increments are shown in Figures III-C-13 and 14. 

4. Plans for Upcoming Period 
This progress report coincides with the completion of 
Nabil A. B. Yehia's Ph.D. studies. At the present time, the 
project is being put on hold until a new Ph.D. candidate with 
an appropriate background for this project can be found. 

Three areas require additional effort to make this 
method a practical tool for micromechanical composite frac- 
ture analysis. They are: 

1) adding iterative steps after each incremental 
solution to control drift from the true solu- 
tion , 

2) developing appropriate criteria to more com- 
pletely represent the fiber-matrix interfaces and 

3) carrying out laboratory experiments to verify 
analytical solutions. 
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Figure III-C-13. Finite Element Mesh for the Second Propagation 
Increment Along the Interface of the Single 
Fiber-Matrix Model. 
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Figure III-C-14. Finite Element Mesh for the Fourth Propagation 
Increment Along the Interface of the Single 
Fiber-Matrix Model. 
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6. Current Publications or Presentations by 
Professor Shephard on this Subject 

"Automated Analysis Model Generation" 

Published in Computer Aided Design in Civil Engi- 
neering, C. N. Kostem and M. S. Sheohard, eds . , 

1984, pp. 92-99. 

"On the Effect of Quarter-Point Element Size", with N. A. B. 
Yehia. 

To be published in the Int. Journal of Num. Meth. 

Engng. 

"Automated Crack Propagation Tracking" 

To be published in Computers and Structures. 

To be presented at Proceedings of Symposium on Ad- 
vances and Trends in Structures and Dynamics, 

Arlington, VA, October 22-25, 1984. 

"Finite Element Model Generation" 

Presented at Siggraph '84, Advanced Topics in - .. 

Solid Modeling, Minneapolis, MN, July 23-27, 1984. 
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III— D DELAMINATION FAILURES OF COMPOSITE LAMINATES 

Senior Investigator: T. L. Sham 

1. Introduction 

Analytical and numerical studies of the delamination 
failures of graphite/epoxy laminates are being investigated. 
The goal of this research is to provide fundamental under- 
standing of, and characterization methodology for, the dela- 
mination failure process in composites. 

2. Status 

A generalized plane strain finite element program has 
been developed for analyzing orthotropic composite layers. 

3. Progress During Report Period 

In order to investigate candidate global parameters that 
could serve as a characterizing parameter for the initiation 
of free-edge delamination, an additional specialized finite 
element code has been incorporated into the finite element 
program. Trial runs on crude finite element meshes have been 
conducted to verify the added code, and production runs on 
meshes based on various lay-ups are about to begin. The par- 
ticular lay-ups are chosen to coincide with those reported in 
the open literature on free-edge delamination experiments. 
Such experimental data would help us to assess the predictive 
capability of the candidate parameters. 
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The development of this finite element program also 
paves the way for investigating innovative approaches in cal 
culating the energy release rate for delamination cracks. 
Typically, the energy release rate is computed by using 
either the virtual crack-closure technique or crack-tip par- 
ameters obtained directly from the singular crack-tip finite 
element. Preliminary analyses indicate that the global 
technique is promising in providing efficient and accurate 
means to calculate the energy release rate. 

Work on applying the "homogenization theory" to compos- 
ite mechanics is also underway, in an effort to clarify the 
effects of higher order constitutive representation of the 
stress and deformation states near free surfaces. 

4 . Plans for Upcoming Period 
The tasks, which are discussed in the preceding section, 
will be continued. 



71 


PART IV 

GENERIC STRUCTURAL ELEMENTS 

IV- A EFFECT OF NOTCH SIZE ON COMPOSITE LAMINATES WITH STRESS 
CONCENTRATIONS 

IV— B IMPROVED BEAM THEORY FOR ANISOTROPIC MATERIALS 
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IV-A EFFECT OF NOTCH SIZE ON COMPOSITE LAMINATES WITH STRESS 
CONCENTRATIONS 

Senior Investigator: D. B. Goetschel 

1. Introduction 

Advanced composite laminates with stress concentrations 
show a notch size effect on failure strength. Strength is 
not governed solely by peak stress but also by the stress in- 
tensity in the vacinity of a peak stress. This was shown 
when tensile specimens with different size holes, which all 

should have the same peak stress concentration factor, failed 

[ 1 ]* 

at significantly different stress levels . This has been 

termed "notch size effect" and can not be predicted by a 

classical stress concentration factor. 

[21 

The "average stress" criterion seems to be the method 
for predicting this effect which is most widely accepted by 
industry. While this method does improve on classical peak 
stress concentration factors, it still has some inaccuracies 
and lacks generality for arbitrary stress fields. This work 
deals with investigating some generalizations of the average 
stress criteria which should make it more versatile, accurate 
and logical. 

Three basic generalizations are being examined. First, 
it is not logical that the average stress calculation be a 
uniform average over a given length without any tapering of 
a weighting function. Therefore, different weighting 

* 

Numbers in brackets in this section refer to the references 
which are listed on page 82. 
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functions will be examined to find a more reliable averaging 
method. Second, the failure analysis will be taken from the 
laminate level down to the ply level to generalize the method 
to irregular geometries. This has already been partially 
verified in previous work. Third, an area average will be 
investigated as opposed to the conventional line average. 

This is more logical, should increase accuracy and help gen- 
eralize the method to irregular geometries and stress fields. 

The proposed modifications to stress concentration fail- 
ure prediction will primarily be investigated experimentally. 
An analytic investigation based on Weibull statistics will 
also be pursued. 

2. Status 

This work was initiated at the beginning of the reporting 
period. Work has concentrated on basic formulation of the 
analysis to be used and on the preparation of the experimental 
program. 

For the analytic formulation, the following three assump- 
tions have been made: 

1) Failure mode is tensile (for now) . 

2) Fiber failure will dominate strength. (Static 
loading only is being considered at this time.) 

3) Notch is free of stress singularities. (It is 
proposed to handle singular elastic stress 
fields by setting a limit on d °x/dx, i.e., 
limit on shear transfer from matrix to fiber.) 
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Finite element analysis will be used to find the stress 
fields for arbitrary geometric configurations . 

3. Progress During Report Period 

During the report period two tasks have been completed. 

An initial finite element computer program for stress 
analysis has been finished in a basic form. Various post- 
processing capabilities still have to be added as the differ- 
ent proposed modifications to the average stress criteria 
are investigated. 

The basic experimental program has been largely com- 
pleted. Published experimental results are not adequate for 
our needs. It was decided that it was important to vary more 
parameters (see below) than has been consistently done so 
far. It was also deemed important to have statistical 
strength information by testing five replicates of each con- 
figuration. Published experimental programs are almost com- 
pletely lacking in this respect. Additionally, most pub- 
lished experimental programs contain incomplete background 
data which makes their utilization difficult for in-depth 
detailed research. 

The experimental program was broken down into five parts 
as follows: 

Part I: We investigated the stacking sequence effects 

( [90/±45/0] , [±45/0/90] and [0/±45/90] ) 
s s s 

and orthotropy effects (two groups: [0^/145] 

a = 1, 2, 4 and [0/±9] , 0 = 15°, 30°, 45°, 60°, 

75°) on failure strength. Although [±45]^ 
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Part II: 


Part III: 


Part IV: 


does not and [90/±30] may not satisfy the 
assumption that fiber failure will dominate 
the failure behavior, it seems worth includ- 
ing both in this testing program because 
[90/±30] s is a "soft" quasi-isotropic lami- 
nate, and [±45 ] 2s is popular for the inves- 
tigation of crack propagation. It is inter- 
esting to have comparisons of [0/90] 2g and 
[±451 ; [ 0/± 45] ; and [0/±60] and [90/ 

£ S S s 

±30] 2s , which might give some indication for 
application of composite structure. Also, 
center cracks and center slots are investi- 
gated by picking up specified orientations 
( [0/±45/90] and [0 2 /±45]); this is for the 
relationship between notch type and material 
properties on failure strength (see Figures 
IV-A-la through lc and Table IV-A-1) . 

We extended the scope. of notch types to in- 
clude double edge crack, double edge semi- 
hole and shoulder (representative for general 
stress raiser cases) [see Figures IV-A-ld 
through If and Table IV-A-2] . 

We also investigated the ability to predict 
strength for more irregular stress distribu- 
tions by cocking cracks and slots at an angle 
(see Figures IV-A-lg and lh and Table IV-A-3) 

Carbide tip twist drilling for holes in com- 
posite laminates is very popular in industry. 
Therefore, the testing data for quasi-iso- 
tropic laminates was compared with corres- 
ponding data in Part I where diamond core 

drills were used. This information will be 
helpful in identifiying how manufacturing 

affects strength (see Table IV-A-4). 




Figure IV-A-1. Specimen Configurations 
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Part V: Strips of different width but equal ratio, — 

w 

(representing the effect of width factor on 
failure strength) were tested (see Table IV-A-5). 

Some remarks can be made based on the experimental test 
results and are listed as follows: 

1) Stacking sequence affects the strength of notched 
and unnotched specimens and in about the same 
proportions . 

2) For center holes with large the average net 
cross-sectional failure stress is about equal to 
the failure stress of unnotched specimens (as 
was expected) . 

3) For quasi-isotropic lay-ups, different types of 
notches with the same nominal dimensions pro- 
duced nearly equal strengths. 

4) Different quasi-isotropic lay-ups, while they 
have the same stiffness behavior, have signifi- 
cantly different strengths. 

5) For the unnotched specimens the Weibull distri- 
bution seemed to be a good model for failure 
probability, with a shape parameter around ten. 

6) No significant specimen width effects were found. 

This indicates that, for widths of at least one 
inch, width is not a critical item and that edge 
effect coupling to stress concentration effects 
is not significant. 

An interesting result from the finite element analysis 
is that the peak stress in the ±45° lay-ups occurs away from 
the notch boundary, in the interior of the specimen. 
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TABLE IV- A- 1 


PART 

ONE OF 

TEST 

PROGRAM: 

CENTER 

HOLES 





Notch Size a 

(in. ) 



Laminate 

0.0 

0.030 

0.065 

0.125 0 

.250 

0.375 0.500 

0.750 

[ 90/±45/0] 

s 

© 1 

NT 11 

© 

iii 

o 

© 

o 

© 

NT 

[±45/0/90] 

s 

© 

NT 

© 

o 

© 

o 

© 

NT 

[0/ ±45/90] ^ 

© 

© 

© 

o 

© 

o 

© 

© 

[0/±45] s 

© 

NT 

© 

o 

© 

o 

© 

NT 

[0 2 /±45] s 

© 

© 

© 

o 

© 

o 

© 

© 

[° 4 /±45] s 

© 

NT 

© 

o 

© 

o 

© 

NT 

[0/±15] s 

© 

NT 

© 

o 

© 

o 

© 

NT 

[0/±30] s 

© 

NT 

© 

o 

© 

o 

© 

NT 

[0/±75] s 

© 

NT 

© 

o 

© 

o 

© 

NT 

[0/±60] 

S 

© 

NT 

© 

o 

© 

o 

© 

NT 

[°/ 90 ] 2 s 

© 

NT 

© 

o 

© 

o 

© 

NT 

[ ±45 1 2s 

© 

NT 

© 

o 

© 

o 

© 

NT 

[90/±30] 2s 

© 

NT 

© 

o 

© 

o 

© 

NT 


1 ©: 5 replicates already tested 

11 'NT: Not tested 

1;L1 o; 5 replicates scheduled for testing 
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TABLE IV- A- 2 

* 

PART TWO OF TEST PROGRAM: ADDITIONAL SYMMETRIC NOTCH TYPES 


Notch Type 

Double Edge Crack 
Double Edge Semi-Hole 
Shoulder 


Notch Size a (in.) 


0.250 0.500 . 0.750 

0 1 © © 

© © © 

© © © 


* 

Laminate : 
1 ©: Tested 


[0/ ±45/90 ] 

[o 2 / ±45 ] s S 


TABLE IV- A- 3 

PART THREE OF TEST PROGRAM: * 

IRREGULAR STRESS FIELDS USING CENTER SLANT CRACKS AND SLOTS 


Slant Angle (6) 
30° 

45° 

60° 

90° 


Laminate: [0/±45/90] s 

[0,/±451 

2 S 

1 @: Tested 


Notch Size a (in.) 


0.250 0.500 0.750 

© X © © 

© © © 

© © © 

© © © 
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TABLE IV- A- 4 

PART FOUR OF TEST PROGRAM: 
MANUFACTURING DRILL TOOL EFFECTS 

Notch Size a (in.) 

Laminate 0 . 030 0.060 0.125 0 . 250 0 . 375 0 .500 

[0/±45/90 ] WT 1 e 11 © © © © 

"HfT: Waiting to be tested 
1;L ©: Center hole drilled by carbide drill 


TABLE IV-A-5 

PART FIVE OF TEST PROGRAM: WIDTH EFFECTS 

Notch Size a (in.) 0.25 0.50 0.75 

Width w (in. ) 0.9 1.8 2.7 

Laminate : © 1 © © 

[0/±45/90] s 


1 ©: Tested 


0.750 

WT 
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4. Plans for Upcoming Period 
The testing program will be completed and the test re- 
sults documented. Some shear test specimens may be included. 
Next, postprocessors for the computer analysis will be writ- 
ten and exercised to test the various proposed generaliza- 
tions to the average stress criterion. Additionally, Weibull 
strength theory will be examined to see if it correlates with 
test results. 


5. References 


1. Waddoups, M. E., J. R. Eisenmann and B. E. Kaminski, 
"Macroscopic Fracture Mechanics of Advanced Composites 
Materials", Journal of Composite Materials, Vol. 5, 1971, 
p. 446. 

2. Whitney, J. M. and R. J. Nuismer, "Stress Fracture Cri- 
teria for Laminated Composites Containing Stress Concen- 
trations", Journal of Composite Materials, Vol. 8, 1974, 
p. 253. 
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IV-B IMPROVED BEAM THEORY FOR ANISOTROPIC MATERIALS 

Senior Investigator: O. Bauchau 

1. Introduction 

Timoshenko beam theory is well-known by structural de- 
signers and is widely used as a first approximation in numer- 
ous structural applications. For solid cross-sectional beams 
made out of isotropic materials, this theory gives accurate 
predictions for aspect ratios L/h >5 (L is the span of the 
beam and h its height). Because it is based on the assump- 
tion that cross sections remain plane after deformation, the 
theory predicts a linear distribution of axial strains. 
However, additional axial strains are induced by warping in- 
compatibilities generated either by specific loading or 
boundary conditions, or by the occurence of nonuniform tor- 
sion or bending. 

The accuracy of Timoshenko beam theory has been shown 
to be rather poor in the case of thin-walled beams because 
of significant warping and shear lag effects and in the case 
of structures made out of highly anisotropic materials. 

The purpose of this work is to investigate the behavior 
of thin-walled box beams made out of anisotropic materials. 

2. Status 

At the beginning of the reporting period the following 
tasks had been completed. 
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a. An analytical solution of the problem was worked out, 
based on the sole assumption that each cross section is in- 
finitely rigid in its own plane but free to warp out-of- 
plane. This assumption has two important implications. 

First, the in-plane displacements of the cross sections are 
fully represented by three rigid body modes, and second, 
any applied transverse load only induces membrane stresses 
in the structure, namely, an axial stress flow, n, and a 
shear stress flow, q. For thin-walled beams, these two 
stress flows are uniform across the thickness of the wall, 
and the other stress components are assumed to be negligible. 

b. A finite element solution of the problem was completed, 
based on the same assumptions outlined above. The results 
for the different approaches were found to be in close 
agreement. 


3. Progress During Report Period 
During the report period two tasks were completed. In 
the first task, the finite element solution was extended to 
allow the treatment of problems presented by structural coup- 
lings. When the axes of orthotropy of the material are not 


aligned with the axis of the beam, such structural couplings 


occur as extension-twisting, extension-bending or bending- 
twisting. 


A specific example is described here to illustrate these 
results. Figure IV-B-1 depicts the thin-walled, rectangular 



40 mm 



Figure IV-B-1. Geometry of the Structure. 



86 


cross-sectional beam that was analyzed. The axis of ortho- 
tropy of the upper laminate was oriented at fifteen degrees 
with respect to the axis of the beam, and the lower laminate's 
was oriented at minus fifteen degrees. The response of the 
structure was first analyzed under a uniform transverse load 
in the y direction. Figures IV-B-2 and 3 show the displace- 
ment, v, in the direction of the load and the twist, 0, 
respectively (note that this twist is due to the coupling 
effect) . For each case, the results of Timoshenko beam 
theory (i.e., neglecting warping effects) were plotted for 
comparison. Due to warping deformations, the tip deflection 
of the beam increased by 12 percent and the tip twist in- 
creased by 78 percent. 

A second case was analyzed where the structure was sub- 
jected to a uniformly distributed torque. Figures IV-B-4 
and 5 show the twist, 0, and the transverse displacement, v, 
respectively (note that this displacement is due to the coup- 
ling effect) . The Timoshenko solution was again plotted for 
comparison. Warping effect considerably increased the tip 
twist of the beam (12%) and the tip deflection (67%) . 

Two conclusions can be drawn from the analyses: a) an 

accurate modeling of the torsional behavior of anisotropic 
beams requires the consideration of warping effects, and b) 
the displacements resulting from structural coupling (i.e., 
twist under transverse load or transverse displacement under 
torque) are drastically influenced by warping deformations. 



Transverse Displacement, v[10 



Figure IV-B-2. Transverse Displacement Along the Span of the Beam 

for hhp Tlni -f- <-\ v-m T /n -• n . . 


for the Uniform Load Case. (P QV . = 10,000 ll/m) 
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Span of Beam 


Figure IV-B-4. Twist Along the Span of the Beam for the Uniform Torsion 

Case. (M = 10,000 N) 
o 
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In the second task, an experimental program was initi- 
ated to assess the accuracy of the theory. The test speci- 
mens were thin-walled box beams with a rectangular cross 
section (see Figure IV-B-6) . The upper and lower skins are 
graphite/epoxy laminates, and the webs are aluminum C-chan- 
nels. The skins are bolted onto the webs to form the box 
beams. The beam was rigidly clamped at one end, while a con- 
centrated transverse load or torque was applied at the other 
end. Two types of specimens were manufactured: First, the 

principal axis of orthotropy of the upper and lower skins 
was aligned with the axis of the beam so that no structural 
couplings were present, and second, the principal axis of 
orthotropy of the upper and lower skins was oriented at an 
angle with respect to the axis of the beam so that structural 
couplings were present. 

The response of the beam was measured by means of the 
following instrumentation. Dial gages were used to measure 
the tip deflection and tip twist of the structure. Strain 
gages were placed near the root of the beam to measure the 
strain distribution across the upper skin. Axial gages, as 
well as rosette gages, were used so that both axial and shear 
strains were measured. 

4. Plans for Upcoming Period 

Effort will be concentrated on the completion of the 
testing program. First, an experimental verification of the 
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warping theory will be sought by comparing measured deflec- 
tions and strain distributions with the predicted ones. 

Another point of interest will be the analysis of struc- 
tural couplings. Various types (i.e., extension-twisting, 
and/or extension-bending and/or bending-twisting) and dif- 
ferent strengths of structural couplings can be obtained by 
varying the configuration of the test specimen. 

Finally, to complete the comparison between analytical 
approaches and experimental results , one specimen configura- 
tion will be modeled using 3-D plate finite elements. 

5. Current Publications or Presentations by 
Professor Bauchau on this Subject 

"A Beam Theory for Anisotropic Materials" 

To be published in the Journal of Applied Mechanics. 
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PART V 


PROCESSING SCIENCE AND TECHNOLOGY 


V-A VARIATION OF RESIN PROPERTIES THROUGH THE THICKNESS OF 
CURED SAMPLES 

V-B NUMERICAL ANALYSIS OF COMPOSITE PROCESSING 
V-C HEAT TREATMENT OF METAL MATRIX COMPOSITES 
V-D INITIAL SAILPLANE PROJECT: THE RP-1 

V-E SECOND SAILPLANE PROJECT: THE RP-2 




97 


V-A VARIATION OF RESIN PROPERTIES THROUGH THE THICKNESS OF 
CURED SAMPLES 

Senior Investigator: B. Wunderlich 

1, Introduction 

A variation of resin properties through the thickness of 
a N,N,N'N' -tetraglycidyl-4 ,4 ' -diaminodiphenylmethane (TGDDM) / 
4 ,4 ' -diaminodiphenylsulfone (DDS) sample, a commonly used 
aerospace epoxy matrix, was initially proposed for study. 

The degree of cross-linking of TGDDM by DDS in a specimen 
was to be determined by quantitatively superposing the glass 
transition region (through the glass transition temperature, 

T , and the change in heat capacity at T , AC (T ) ) of an 

y y . p y 

uncharacterized material with that of a characterized refer- 
ence. Different specimens of the same sample were to be 
analyzed to assess heterogeneity. 

Curing of TGDDM/DDS systems involves two types of pos- 
sible cross-linking reactions: epoxide-epoxide and epoxide- 

amine. Both types of cross-linking reactions proceed via a 

[i] * 

ring-opening reaction of the epoxide . For superposition 

to be possible, the epoxide-amine reaction should primarily 

determine the extent of cross-linking, in order that T be 

y 

dependent upon the amount of hardener (DDS). Structures, 
descriptions and reactions of the components of the TGDDM/DDS 
system are given in [2] . 


Numbers in brackets in this section refer to the references 
which are listed on page 108. 
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2. Status 

Epoxy resins containing various amounts of hardener 

(DDS) were cured in a muffle furnace at 473°K for seven hours. 

[3] 

The TGDDM manufacturer recommends this cure time . T g ver- 
sus ratio-percent of hardener to epoxy resin was measured 
(see Figure V-A-l) . A limit of approximately 475°K was ra- 
pidly reached at only two percent hardener, a 209°K increase 
over the uncured pure resin. Pure TGDDM was cured in a muf- 
fle furnace under the same conditions as the resin/hardener 
system. A T g of 386°K was obtained, an increase of 121°K 
over the uncured pure resin or, approximately, half the in- 
crease seen for the samples with more than two percent hard- 
ener. In addition, Fourier Transform Infrared Spectroscopy 
(FT-IR) data evinces the appearance of ether groups for 
TGDDM/DDS systems cured at 450°K for 2.5 hours for 0 to 35 
weight-percent DDS^] Thus, determining the extent of amine 
cross-linking for the TGDDM/DDS system by superposition of 
the glass transition temperature of a characterized reference 
sample under ordinary cure conditions is complicated due to a 
masking effect of the epoxide-epoxide curing reaction. Pos- 
sible superpositon of AC (T ) may allow differentiation be- 

ir y 

tween the different cures for the TGDDM/DDS system. 

The glass transition region separates glassy and liquid 
behavior for amorphous materials. Therefore, the knowledge 
of the behavior and characteristics of the temperature re- 
gions below, higher and through which the glass transition 



< K ) 
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occurs is critical, in the determination of application par- 
ameters for amorphous materials. Thermal properties within 
the glass transition region; the glass transition temperature, 
T ? the change in heat capacity at T , AC (T ) ; the heat ca- 
pacity hysteresis, AH h ; the broadness of the inner and outer 
glass transition regions, A^ and AT 2 and, additionally, the 

heat capacity of the glass and melt have been studied for 

[41 

amorphous noncross-linked materials . For amorphous, cross- 

linked materials few of such studies exist. 

Initial work by Ueberreiter and Kanig in 1947 using 

volume -temperature measurements on cross-linked styrene- 

divinyl benzene (S-DVB) demonstrated that T g increased with 

[51 

increasing cross-linking . Further analysis of their data 
by Gians and Turner evinced a linear relationship between T g 
and the cross-link density, except at low concentrations of 
the cross-linking agent Additionally, Gians 'and Turner 
found that AH h decreases with increasing cross-linking 1 . 
However, very little further investigations of the interre- 
lationships between T g , AH h and AT (AT 2 is more difficult to 
determine) with increasing cross-linking will be discussed 
in a report which is being prepared for publication. It is 
intended to provide the first step in a fuller understanding 
of thermal behavior as a function of cross-linking. 



101 


3. Progress During Report Period 

a. Systems Investigated 

As of this date we have investigated three cross-linked 
systems: 

1) polysep, commercial cross-linked polystyrenes, 

2) styrene, radically bulk polymerized, with di- 
vinyl benzene as cross-linking agent and 

3) N, N,N' ,N'-tetraglycidyl-4 ,4 "-diaminodiphenyl- 
methane ( TGDDM) . 

Polysep systems are commercial (Polysciences) S-DVB co- 
polymer beads containing between 0.5 and 12 percent-by-weight 
DVB. The second type of S-DVB copolymer used was synthesized 
by bulk polymerization with 0.5 percent-by-weight of benzoyl 
peroxide. Finally, TGDDM, a tetraf unction epoxy, was ther- 
mally cured at 473°K, by itself, via a ring opening reaction. 

b. Glass Transition Temperature 

All systems were found to show an increase in T g with 
an increase in cross-linking except for S-DVB samples at low 
concentrations of DVB (see Figures V-A-l, 2, 3). To investi- 
gate this trend at low concentrations of DVB further, various 
samples of bulk polymerized S-DVB, containing low concentra- 
tions of DVB, were made. A minimum two percent-by-weight of 
DVB was noted for both types of S-DVB samples investigated, 
and a lower maximum at 0.25 percent of DVB. Dissolution ex- 
periments in toluene demonstrated that the 0.25%-DVB sample 
was soluble, while samples with increased DVB content formed 
increasing amounts of insoluble gels. 




LEGEND 


Literature data for a heating 
rate of 20°K/min [6] 

Polysep samples 

Literature data for a linearly 
extrapolated zero heating rate 

Bulk polymerized samples 


i i i 

5 30 35 40 45 50 

i DVB 

ures as a Function of DVB Concentration 



T a /K 


405 



Wt.% DVB 

Figure V-A-3. Low Concentration Region of Figure V-A-2 
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The behavior of the S-DVB samples studied may be ex- 
plained in light of a heterogeneous structure. In general, 

many cross-linked materials are thought to consist of a het- 

[7] 

erogeneous network structure . Electron microscopy studies 

of cross-linked materials have revealed nodular-like or 

colloidal-like structures of high cross-link density embedded 

[8 9 101 

in a matrix of lesser cross-linking ' ' i Solomen et 
ai^ll] j^ ave SU gg es ted that these "nodules" are formed early 
in the curing process, competing with the pure polymerization 
process, and then they grow. Competition between the forma- 
tion of dispersed cross-link points in the matrix, nodule 
formation and growth may give rise to the "roller coaster" 
like relationship seen in Figure V-A-2 . That is, highly dis- 
persed cross-links may hold the matrix more rigidly than the 
more widely space "nodules" of approximately the same DVB 
content. 

c. Heat Capacity Hysteresis Effect 

The decrease of hysteresis with increasing cross-link 
density was investigated using S-DVB samples as a model sys- 
tem to obtain quantitative information on the structure of 
cross-linked materials. A relationship between the broad- 
ness of the glass transition region, AT X , and the disappear- 
ance of the hysteresis peak, AH h , was noted (see Figure V- 
A-4) . Such relationships have been observed in block copol- 
ymers and blends, although reduction of the hysteresis has 
also been noted in semicrystalline and mesophase 




AH h /J mol" 



Figure V-A-4. Change of the Glass Transition Hysteresis Peak Area, AH h , and the Breadth of the Glass 
Transition Region, AT^, as a Function of Weight-Percent DVB for Cross-Linked Polystyrene. 
(Cooling rate 0.5 K°/min ,. heating rate 20°K/min, Tg determined by the ^ACp meth- 
od from a 20°K/min cooling and heating) 
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materials ' '*' 3 A similar behavior was also found 

for TGDDM (see Figure V-A-5) . 

In general, with increasing cross-linking, T g was found 
to increase, the glass transition to broaden and the hyster- 
esis peak to be reduced. Exceptions to the previous gener- 
ality occur at low concentrations of DVB (less than 2% by 
weight). Increasing feeds of DVB cause a decrease in T g , and 
with this decrease, the glass transiton narrows and the hys- 
teresis increases (see Figures V-A-2and'4). Further feeds 
of DVB result in a reversion to the original trend. At very 
high degrees of cross-linking, TGDDM samples also show an ex- 
ception to the previous generality. That is, with increasing 
cross-linking at high cross-link densities, even though T g 
still increases, the glass transiton region begins to narrow, 
and the hysteresis peak increases again (see Figure V-A-5) . 

d. Broadening of the Glass Transition 

The regions on both sides of the glass transition tem- 
perature were compared to assess the broadening present in 
cross-linked samples. Different comparisons were made on the 
T -T and T -T regions and on the T -T and T -T regions. 

IB E 2 g B Eg 

The type of broadening seen for both the S-DVB and TGDDM sam- 
ples was symmetrical. That is, for the same thermal history, 
the differences between T E ~T 2 and T^-Tg or T E ~T g an( ^ T g -Tg 
were approximately the same for each of the samples studied 
at various degrees of cross-linking. 




Figure V-A-5. Change of the Glass Transition Hysteresis Peak Area, AH h , and the Breadth of the Glass 
Transition Region, AT^, as a Function of Glass Transition Temperature T , for 
TGDDM Cured by Itself. (Cooling rate 0.5°K/min, heat rate 20°K/min, T g g de- 
termined by the 4AC method frpm a 20°K/min cooling and heating) 
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Generally, the increasing breadth of the glass transi- 
tion with increasing cross-linking may be due to the hetero- 
geneous nature of the samples examined. Differing distances 
between cross-link points and "nodules" should give rise to 
inhomogeneities of differing rigidities. Since a more rigid- 
ly held polymeric chain has been related to increasing the 
T one would expect, as the amount of cross-linking is 

increased, both the T g and the distribution of inhomogenei- 
ties would increase and, hence, broaden the glass transition. 

4. Plans for Upcoming Period 
Heat capacity data and AC (T )'s for cross-linked ma- 

ir y 

terials will be collected, evaluated and presented in our 
next report. Comparisons to noncross-linked materials will 
be made. A full interpretation of the many observed experi- 
mental facts is planned in parallel with the completion of 
the thesis of Mr. L. Judovits. 
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V-B NUMERICAL ANALYSIS OF COMPOSITE PROCESSING 

Senior Investigator: II. S. .Shephard 

1. Introduction 

The complexities in both design and manufacture of com- 
posite materials demand the use of electronic computation to 
efficiently utilize these materials. Presently, the analyses 
performed in connection with the production of a new part are 
related to proper design for use in a service environment. 

In the future, there will be a tendency to perform more anal- 
yses associated with the manufacture of a product. The goal 
of this new project is to develop the appropriate analysis 
tools needed to describe the processing of continuous fiber 
resin matrix composites. 

Composite properties are directly related to the quality 
of the lay-up that results from the fabrication process. 

This manufacturing process is a very complex operation that 
depends on many variables. An understanding of the fabrica- 
tion processes will allow proper control of the manufacturing 
variables to obtain desirable qualities in the finished pro- 
duct. Numerical analysis provides a means for the solution 
of the complex governing differential equations that describe 
composite processing. The ability to solve the governing 
equations, coupled with the proper input data and comparison 
with experimental results, will lead to a better understanding 
of composite processing. 
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2 . Status 

The systems that are of interest for this initial work 
include both thermosetting and thermoplastic resins. Pro- 
cessing will take a slightly different general form for each 
of these types of matrices. However, there is a good deal 
of overlap in the types of numerical analyses that are re- 
quired to understand the curing of thermoset composites and 
the processing of thermoplastic composites. For these anal- 
yses, there is a considerable amount of input data that is 
required, and at this point, it is assumed that the proper 
information will be available. A description of the pro- 
cessing and analyses of both thermosets and thermoplastics 
follows. 

a. Thermosets 

[ 1-41 * 

Springer and associates published several articles 

dealing with the analysis of the curing process for epoxy 
matrix composites. In their work, the curing process was ex- 
amined with a series of one dimensional models for the study 
of heat generation and transfer, resin flow, void volume 
changes and thermal stresses. These analyses were carried 
out on a macromechanical level and are limited to a re- 
stricted geometry (infinite plate) . Their work represents a 
good first step in the analysis of the curing process but is 
lacking in a few areas. However, the basic approach is valid 

and is described in the following paragraphs. 

* 

Numbers in brackets in this section refer to the references 
which are listed on page 120. 
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There are several coupled processes that occur during 
curing that must be modelled. First, there is a thermal- 
chemical process that deals with heat transfer and the heat 
generated from the exotherm. This process is coupled to the 
flow process by convection heat transfer and by flow during 
the reaction or just flow with a nonuniform temperature dis- 
tribution. The flow model is concerned with the flow of a 
viscous fluid through an orthotropically porous medium. In 
addition, a void model is necessary to account for the vola- 
tiles produced by the exotherm and the voids that are mechan- 
ically introduced. The void model is coupled to the thermo- 
chemical model by the temperature and the degree of completion 
of the reactions. Finally, there is a thermo-mechanical 
model to describe the thermal stress distributions that arise . 

The thermal-chemical model is a heat transfer model 
that has terms relating to heat storage, heat generation and 
heat transfer. The heat stored in the material is related 
to the density and the heat capacity of the material. The 
heat generated is related to the degree of completion of the 
cross-linking chemical reactions. The heat transfer term 
accounts for heat flowing by means of conduction, convection 
and radiation. 

The flow model requires data for the viscosity of the 
resin through the curing process. This model must account 
for flow of the resin in all directions both parallel and per- 
pendicular to the fibers. The analysis of the flow model will 
result in knowledge of the resin content of the final lay-up. 
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The void model must account for the generation of gases 
during the exotherm. The volume of trapped gas has to be re- 
lated to the amount of gas present, the viscosity of the 
resin and the pressure applied to the lay-up. Analysis of 
the void model will indicate the volume fraction of voids 
that results from the curing process. 

The thermal-mechanical model relates the temperature 
and degree of cure to the residual thermal stresses. The 
model must account for the change in the resin's properties 
from those of a viscous fluid to those of a viscoelastic 
solid. This analysis must track the stresses through the 
change from viscous fluid to solid. 

b. Thermoplastics 

Similar analyses for thermoplastic composites have not 
been carried out. The processing of thermoplastics is more 
difficult because of the non-Newtonian behavior of the resin 
and the high temperatures and pressures required to obtain a 
viable engineering material. There are several methods for 
processing thermoplastic composites which are all character- 
ized by the flow of the matrix material, or a solution of the 
matrix material, into dry fibers. 

One method for forming thermoplastic composites begins 

by laying up alternate layers of thermoplastic film and gra- 
[ 5 ] 

phite cloth . Heat and pressure are then applied to complete 
the forming process. The lay-up is heated above the glass 
transition temperature of the resin in order to reduce the 
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viscosity. Pressure is then applied to generate the flow of 
the resin into dry fibers. A second processing step might 
be carried out from this point to form the thermoplastic com- 
posite into the desired shape with heat and pressure. 

The thermal processing for thermoplastics is basically 
different from the curing process of a thermoset due to the 
absence of a chemical reaction. However, the energy equation 
includes a term for the heat of crystallization which is 
treated in the same manner during numerical analysis as the 
heat of reaction term. In addition, there is a possiblity 
for thermal expansion anisotropy in the resin and shear 
heating of the resin. 

The flow process for thermoplastics is more complicated 
than the parallel process for thermosets for two reasons. 

The first is linked to thermal stresses caused by differences 
in thermal coefficients of expansion of the resin and fibers. 
These stresses must be tracked through the transition of the 
resin from a viscoelastic fluid to a viscoelastic solid. A 
second source of residual stresses is that associated with 
the dimensional change upon crystallization of the matrix. 

3. Progress During Report Period 

To obtain a feel for the type of numerical analyses that 
are going to be required for this project, a one-dimensional 
analysis code was implemented to solve the parabolic energy 
equation associated with the curing of a thermosetting 
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material. This analysis utilizes finite elements for the 
spatial domain and finite differences for the time domain. 

The general macroscopic energy equation, neglecting 
kinetic energy, that is valid for both thermosets and thermo 
plastics is given by: 


PD(CT) 

Dt 


V* [K] VT + Q(T,t) 


where: p = the composites density 

t = time 
T = temperature 

C = the composites heat capacity 

Q = rate of heat energy due to internal 
sources or sinks 

[K] = thermal conductivity matrix 


and for cartesian coordinates 
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The Q term is used to represent the heat of reaction, heat 
of crystallization or viscous heating. 

Several assumptions are made in the solution of the en- 
ergy equation for the curing of a thermoset composite; the 
composite being cured is a thin infinite plate, convected 
heat transfer terms are negligible, and there are no viscous 
losses. For this initial analysis, the coefficients in the 
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equation are taken to be constant for position, time and 
temperature. With these assumptions, the energy equation 
reduces to the nonlinear parabolic equation given by: 

pc|f= k— + Q(T,t) 

The analysis code was written using linear finite ele- 
ments to discretize the thickness of the laminate. A forward 
difference scheme is used for time integration. This scheme 
was chosen due to its simplicity in handling the nonlinearity 
associated with the heat of reaction term and its ease of 
implementation. The program accepts Dirichlet boundary con- 
ditions that are functions of time. 

The program was tested in two ways. First, with an ex- 
ample that did not contain any internal heat sources and thus 
had an exact analytic solution, and second, with an example 
that included an expression for the heat of reaction and did 
not have an exact solution. For the purposes of testing, the 
heat of reaction term was represented by two uncoupled Arrhen- 
ius rate equations. Two equations were used to represent the 

two major reactions that occur when curing Hercules 3501-6 

[ 2 ] 

resin . The coefficients of the rate equations were chosen 
such that good agreement was obtained with the data of Refer- 
ence [2] . The boundary conditions used for this problem are 
the same as those used in References [3] and [4] and consist 
of a ramp temperature increase from room temperature to 177 °C 
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in 60 minutes. The analysis for this problem yielded solu- 
tions for the temperature and degree of completion of the 
reactions as functions of position and time. Figure V-B-l 
shows the temperature as a function of the thickness at a 
particular instant of time. The temperature in the center of 
the laminate is higher than that at the boundaries due to the 
exotherm. This example was used to demonstrate that the code 
was working properly and not necessarily to depict reality. 

In order to obtain realistic results, a more accurate descrip- 
tion of the kinetics of the reaction is required. Inclusion 
of an improved reaction model will only require minor code 
modifications . 


4. Plans for Upcoming Period 
The one dimensional solution of the energy equation for 
thermosets will be expanded to incorporate all sources of 
nonlinearity. The coefficients will be allowed to vary with 
position in the laminate, degree of cure and temperature. 

The literature will be examined for experimental results that 
can be used to test this model and determine the importance 
of these terms. 

The finite element code ABAQUS can be used to test the 
need for analysis of the energy equation with constant coef- 
ficients in more than one dimension. ABAQUS has solid and 
shell heat transfer elements that can handle the heat of re- 
action with user supplied subroutines. 
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Depth 


Figure V-B-l. Temperature Profile Through Thickness. 



120 


Solution of the flow equation will begin by finding the 
constitutive relations necessary to capture the important 
features of the non-Newtonian behavior of the candidate poly- 
meric matrices. 
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V-C HEAT TREATMENT OF METAL MATRIX COMPOSITES 

Senior Investigator: N. S. Stoloff 

1. Introduction 

During this report period a new program on metal matrix 
composites produced by directional solidification was begun. 
The objective of this work is to utilize postsolidification 
heat treatments to improve the mechanical properties of car- 
bide-reinforced eutectic composites. The alloy system chosen 
for study is Fe-Mn-Cr-C, which can be directionally solidi- 
fied at unusually high rates (to 25 cm/hr) and still maintain 
its composite structure. A further advantage of this system 
is that in the as-solidified condition it is one of the 
strongest iron-base alloys known; its strength at elevated 

temperatures is comparable to that of iron-nickel and cobalt- 

* 

base superalloys (see Figure IV-C-1) . As such it is a face 
centered cubic Fe-Mn-Cr solid solution, and the fibers are 

Cr 7 C 3 * 

The specific properties to be measured as a function of 
postsolidification heat treatment are tensile strength and 
fatigue resistance, in both longitudinal and transverse ori- 
entations . 


Lemkey, F. D. , E. R. Thompson, J. C. Schuster and H. Nowotny, 
"In Situ Composites IV", F. D. Lemkey, H. E. Cline, M. 

McLean, eds.. North Holland, NY, 1982, p. 31. 
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Figure V-C-l. Fe and Co Alloys Designed for High Temperature Usage 
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2. Progress During Report Period 
Immediately upon commencing this work, a directionally 
solidified ingot of Fe-20%Cr-10Mn-3.4C was obtained from Dr. 
Frank Lemkey of United Technologies Research Center. This 
material had been solidified at 25 cm/hr. 

Exploratory heat treatments to develop increased strength 
and hardness were carried out. All treatments consisted of a 
high- temperature solution anneal , followed by a low tempera- 
ture age. Hardness versus aging time was measured for each 
heat treatment. Optimum hardness was obtained by solution- 
ing the alloy at 1170C 0 , followed by again at 850 °C for eight 
hours. The results of that experiment are shown in Figure 
V-C-2. The microstructures before and after aging are shown 
in Figures V-C-3a and 3b, respectively. Aging appears to 
have developed relatively fine carbides dispersed in the ma- 
trix between the massive carbide fibers. 

Subsequently, small compression specimens, 3/16-inch 
diameter by 3/8-inch height were machined from the as-DS rod; 
half of these were heat treated to peak hardness. Initial 
testing of the as-DS material revealed a yield strength near 
270 ksi. The load required to deform this sample was so high 
that new samples of smaller diameter need to be machined in 
order to measure tensile properties of aged material. 

A second alloy, Fe-30Cr-3C, solified at 25 cm/hr, has 
recently been received from Dr. Lemkey. 




Figure V-C-2. Hardness as a Function of Aging Time for the Test Specimen. 


124 




a) As-Received, x 500 



Figure V-C-3. Microstructure of 



b) 1170 P C for 6h, Aged 85Q P C for 24h, 

-3.4C x 1000 


Specimen Before and After Aging. 
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3. Plans for Upcoming Period 
New compression samples are to be machined so that tests 
at both 25°C and at temperatures to 850°C can be conducted. 
These tests will serve as the basis for later tensile and 
fatigue studies on larger ingots of this alloy. 

The second ingot, Fe-30Cr-3C, will be examined for op- 
timum heat treatment response and then prepared for compres- 
sion testing as described above. 


4. Current Publications or Presentations by 
Professor Stoloff on this Subject 

"Eutectic Composites for High Temperature Applications" 

To be presented at the Joint U.S. /Japan Seminar 
on Superalloys, Mt. Fuji, Japan, December 7-11, 

1984. 

"Current Status and Prospects of Eutectic Composite Super- 
alloys" 

To be presented at the Superalloy Seminar, Chung 
Shan Institute of Science and Technology, Lung-Tan, 
Taiwan, December 17-18, 1984, 
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V-D INITIAL SAILPLANE PROJECT: THE RP-1 

Senior Investigators: F. P. Bundy 

R. J. Diefendorf 
H. Hagerup 

During this reporting period the RP-1 glider has been 
stored on the balcony of the Composite Materials Shop Area 
in the Jonsson Engineering Center under ambient conditions 
of temperature and humidity. 

The glider was assembled inverted and given its "annual" 
static strength test on November 15, 1984 as a simple bending 
wing loading up to over 4 G's, in the same manner as in pre- 
vious years. The deflections at the wing tips due to the 
staged wing loadings are compared, for the four years of 
testing, in Table V-D-l. 


TABLE V-D-l 

ANNUAL STATIC STRENGTH TESTS FOR THE RP-1 


Wing Tip Deflection (cm) 


buau 





(lbs/wing) 

1984 

1983 

1982 

1981 

0 

0 

0 

0 

0 

150 

5.72 

5.48 

5. 8 

5.7 

300 

11.45 

10.87 

11.7 

11.3 

430 

15.17 

14.50 

15.4 

15.5 

550 

18.15 

17.26 

18.4 

18.1 

590 

18.56 

17.72 

— 

— 

550 

18.22 

17.40 



430 

15.49 

14.72 



300 

11.83 

11.20 



150 

6.18 

5.80 



0 

.40 

.25 
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Although an attempt was made to make the test setup the 
same for each year, there were some minor differences in the 
supporting system which could account for the variations 
noted in the table for the same wing loadings. The overall 
conclusion is that the airframe strength and rigidity are 
not degrading with time, and that the structure is airworthy. 
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V-E SECOND SAILPLANE PROJECT: THE RP-2 


Senior Investigators: F. P. Bundy 

R. J. Diefendorf 
H. Hagerup 


1. Status 

At the last report period the air frame with the re- 
built wings and redesigned and rebuilt wing beam linkage had 
been tested to over +6 G's in simple bending and to a nega- 
tive torque of 7000 pounds with completely acceptable results, 
except for some forward bending of the wing which caused 
about a 3/16-inch pullout of the aft wing pins from the soc- 
kets in the fuselage. Also, it was reported that a design 
had been developed which would provide tension linkage be- 
tween the aft pin areas of the wings and prevent pullout 
of the aft pins during heavy loading, and that a full size 
model of this tension linkage had been constructed and force- 
tested. 

2. Progress During Report Period 
During the current report period the aft wing tension 
linkage system was built into the wings and fuselage. After 
the installation, the air frame was assembled inverted and 
put through the regular +6-G sand bag loading test, as be- 
fore. This was the third +6-G loading test to which the re- 
built wings and wing beam linkage had been subjected. It 
behaved very satisfactorily, including negligible pull out 
of the aft wing pins. On the basis of the static strength 
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and deflection tests, the wing/fuselage assembly is capable 
of any aerodynamic loading within conditions of the flight 
envelope for the aircraft. 

The towhook/fuselage structure was force-tested to the 
FAA standards of 2 G's in directions ranging from 30-degrees 
up, 75- degrees down and 30-degrees sidewards. The RP-2 
structure held up satisfactorily, and the towhook release 
worked under those specified conditions of loading. 

The FAA-specified static strength tests of the empen- 
nage/fuselage structure were done and passed satisfactorily; 
likewise for the landing gear/fuselage structure. 

After all the above-mentioned major strength tests were 
accomplished, work proceeded on fabrication of the foredeck 
and rudder pedal assembly, fitting and mounting arrangements 
for the fuselage skins and a control stick trim lever assembly. 

3. Plans for Upcoming Period 
During the next two reporting periods it appears that 
the glider can be completed and made ready for flight testing. 
The work remaining is mainly the finishing of the fuselage 
shell, cockpit rim, canopy rim frame instrument panel and 
painting. 

4. Current Publications or Presentations by 
Professor Bundy on this Subject 

"Composite Materials, Engineering and Gliders at RPI" 

Published in 1983-1984 Hangar Flying, by invita- 
tion of the editor of the Flying Engineers 


PART VI 

TECHNICAL INTERCHANGE 
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Technical meetings, both on- and off-campus, provide 
for the interchange of technical information. In order to 
assure that a Rensselaer faculty/staff member can participate, 
a central listing of upcoming meetings is compiled, maintained 
and distributed on a quarterly basis. The calender for this 
reporting period is shown in Table VI-1. Table VI-2 shows 
the meetings attended by RPI composites program faculty/ staff/ 
students during the reporting period. Some on-campus meet- 
ings, with special speakers particularly relevant to compos- 
ites, are listed in Table VI-3. A list of composites-related 
visits to relevant organizations, attended by RPI faculty/ 
staff/students, along with the purpose of each visit is pre- 
sented in Table VI-4. 

The diversity of the research conducted within this pro- 
gram has increased over the last several years. To insure 
information transfer, once-a-week luncheon programs have been 
held among the faculty and graduate students involved (listed 
in Part VII. Personnel - of this report). These meetings 
are held during the academic year and are known as "Brown Bag 
Lunches" (BBL's), since attendees bring their own. Each BBL 
allows an opportunity for graduate students and faculty to 
briefly present plans for, problems encountered in and recent 
results from their individual projects. These seminars also 
are occasions for short reports on the content of off-campus 
meetings attended by any of the faculty/staff participants 
(see Tables VI-2 and VI-4) and for brief administrative 
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TABLE VI- 1 

CALENDAR OF COMPOSITES -RELATED MEETINGS 
(April 30, 1984 through September 30, 1984) 


5/1-4 

5/14-16 

5/16-18 

5/29-30 

6/10-15 

6/25-29 

7/2-6 

7/4-5 

7/8-14 

7/11-13 

7/16-20 

7/19-20 

7/23-27 

7/28-8/4 

8/2 

8/7-9 

8/13-16 


AIAA Annual Meeting and Technology Display, Washing- 
ton, DC. Sponsored by AIAA. 

25th Structures, Structural Dynamics and Materials 
Conference, Palm Springs, CA. Sponsored by AIAA/ 
ASME/SAE . 

AHS Annual Forum and Technology Display, Washing- 
ton, DC. Sponsored by AHS. 

USAF Aeronautical Systems Division Advisory Group 
Meeting, Dayton, OH. Sponsored by USAF. 

5th International Conference on Experimental Mechan- 
ics, Montreal, Canada. Sponsored by SESA/BSSM. 

Composites and Paper Physics Technology Interchange, 
Stockholm, Sweden. Sponsored by NSF. 

Carbon 84 Conference, Bordeaux, France. Sponsored 
by the University of Bordeaux. 

Conference on Materials in Engineering, Solihull, 
England. Sponsored by PRI. 

1st Gordon Conference on Liquid Crystals, New 
London, NH. Sponsored by Gordon Conferences. 

JSPE International Symposium on Design and Synthesis, 
Tokyo, Japan. Sponsored by JSPE. 

10th International Liquid Crystal Conference, York, 
England. Sponsored by the University of York. 

Materials Research Council Meeting, La Jolla, CA. 
Sponsored by DARPA and MRC. 

Siggraph '84, Advanced Topics in Solid Modeling, 
Minneapolis, MN. Sponsored by Siggraph. 

Experimental Aircraft Association Annual Meeting, 
Oshkosh, WI. Sponsored by EAA. 

IUPAC Working Party on Polymer Characterization and 
Rheology Meeting, Wilmington, DE. Sponsored by IUPAC. 

17th ASTM National Symposium on Fracture Mechanics, 
Albany, NY. Sponsored by ASTM. 

ACEE Composite Structures Technical Conference, 
Seattle, WA. Sponsored by ACEE. 
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TABLE IV- 1 continued 


8/13-17 

8/19-25 

8 / 20-22 

8 / 20-22 

8/21-23 

8/26-31 

8/29-31 

9/2 

9/4-7 

9/9-11 

9/17-21 

9/17-21 

9/17-21 

9/18-20 

9/23-26 

9/24-28 


Army Sagamore Conference, Lake Lucerne, NY. Spon- 
sored by the United States Army. 

16th International Congress of Theoretical and Ap- 
plied Mechanics, Copenhagen, Denmark. Sponsored by 
AIAA/IUTAM. 

Experimental Mechanics of Fiber Reinforced Composite 
Materials, Dayton, OH. Sponsored by SESA. 

Atmospheric Flight Mechanics Conference, Seattle, 

WA. Sponsored by AIAA. 

Triservice Workshop on Dynamics, Mechanical Proper- 
ties, Characteristics and Processes of Lightweight 
Metal Matrix Composites, Leesburg, VA. Sponsored 
by ARO/ONR. 

ACS 188th National Meeting, Philadelphia, PA. 
Sponsored by ACS. 

Colloqium on Mechanical Characteristics of Load- 
Bearing Fiber Composite Laminates, Brussels, Bel- 
gium. Sponsored by EUROMECH. 

IUPAC Working Party on Polymer Characterization and 
Rheology Meeting, Wilmington, DE. Sponsored by IUPAC. 

Advanced Research Workshop on Applications of Frac- 
ture Mechanics to Cementitous Composites, Evanston, 

IL. Sponsored by Northwestern/ARO . 

2nd Workshop on Advanced Materials : Recent Advances 

in Nondestructuve Evaluation of Materials, Traverse 
City, MI. Sponsored by AIAA. 

Colloqium on Mathematical Programming Methods for 
the Plastic Analysis of Structures, London, England. 
Sponsored by EUROMECH. 

Composite Materials: Processing, Quality Assurance 

and Repair, Los Angles, CA. Sponsored by UCLA. 

Failure Modes in Composites VII: Deformation Asso- 

ciated with Fracture, Fall Meeting, Warrendale , PA. 
Sponsored by TMS/AIME. 

International Conference on Specialty Polymers: 

Present and Future, Birmingham, England. Sponsored 
by BS . 

13th NATAS Meeting, Philadelphia, PA. Sponsored by 
NATAS . 

International Symposium in Composites: Materials 

and Engineering, Newark, DE. Sponsored by the Uni- 
versity of Delaware. 
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TABLE VI-2 

COMPOSITES -RELATED TECHNICAL MEETINGS ATTENDED OFF-CAMPUS 
(April 30, 1984 through September 30, 1984) 


5/16-18 

6/25-29 


7/2-6 


7/8-14 


7/16-20 


7/19-20 


7/23-27 


AHS Annual Forum and Technology Display (Prof. Die- 
fendorf) , Washington, DC. 

Composites and Paper Physics Technology Interchange 
(Prof. Sternstein) , Stockholm, Sweden. 

Professor Sternstein presented the invited 
lecture: "Mechanical Characterization of 

Composites " . 

Carbon '84 Conference (Prof. Diefendorf) , Bordeaux, 
France. 

Professor Diefendorf presented the invited 
lecture: "Chemical Vapor Disposition of 

Carbon" and the papers: "Characterization 

of Isotropic Phase and Mesophase Pitch Frac- 
tions", Effects of Lithium Radiation Tech- 
niques on the Determination of Pitch Molec- 
ular Weight" and "Pitch Solvent Interactions 
and Their Effects on Mesophase Formation" . 

1st Gordon Conference on Liquid Crystals (Prof. 
Wunderlich) , New London, NH. 

Professor Wunderlich presented the invited 
lecture: "Condis Crystals" and the paper: 

"Condis Crystals of Linear Macromolecules". 

10th International Liquid Crystal Conference (Prof. 
Diefendorf) , York, England. 

Professor Diefendorf presented the paper: 
"Mesophase Formation in Polynuclear Aromatic 
Compounds " . 

Materials Research Council Meeting (Prof. Stern- 
stein) , La Jolla, CA. 

Professor Sternstein presented the invited 
lecture: "Mechanical Characterization and 

Failure of Composites". 

Siggraph '84, Advanced Topics in Solid Modeling 
(Prof. Shephard), Minneapolis, MN. 

Professor Shephard presented the paper: 

"Finite Element Model Generation". 
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TABLE VI - 
7/28-8/4 

8/13-16 

8/13-17 

8/26-31 

9/2 

9/23-26 

9/24-28 


2 continued 


Experimental Aircraft Association Annual Meeting 
(Prof. Diefendorf) , Oshkosh, WI. 

Professor Diefendorf presented the paper: 
"Fabrication of Composite Materials at RPI". 

ACEE Composite Structures Technical Conference 
(Prof. Loewy) , Seattle, WA. 

Army Sagamore Conference (Prof. Diefendorf) , Lake 
Lucerne, NY. 

Professor Diefendorf presented the paper: 
"Materials Characterization of Carbon Fi- 
bers" . 

ACS 188th National Meeting (Prof. Sternstein) , 
Philadelphia, PA. 

Professor Sternstein presented the paper: 
"Viscoelastic Characterization of Compos- 
ites" . 

IUPAC Working Party on Polymer Characterization and 
Rheology Meeting (Prof. Sternstein), Wilmington, DE. 

13th NATAS Meeting (Prof. Wunderlich) , Philadel- 
phia, PA. 

Professor Wunderlich gave the lecture: 

"The Glass Transition of Poly (ethylene 
terephthalate) " . 

International Symposium in Composites: Materials 

and Engineering (Prof. Diefendorf), Newark, DE . 
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TABLE VI -3 

COMPOSITES-RELATED MEETINGS/TALKS HELD AT RPI 
(April 30, 1984 through September 30, 1984) 


Topic 

"Geometric Modeling for 
Mechanical Analysis and 
Synthesis" 

"Structural Optimization: 
Toward a Reliable CAD Tool" 


"Coalescence of Surface 
Fatigue Cracks" 


"Ductile Fracture" 


"An implicit Time Integra- 
tion Scheme for Inelastic 
Constitutive Equations with 
Internal State Variables" 


Date Speaker (s ) 

5/1 Dr. M. Y. Zarrugh 

University of Michigan 


Dr. C. Fleury and 

5/31 Mr. V. Braibant 

University of Liege 
Belgium 

7/23 Prof. K. Iida 

University of Tokyo 
Japan 

8/8 Dr. J. W. Hancock 

University of Glasgow 
Scotland 

9/11 Prof. F. G. Kollman 
Technishe Hochschule 
Darmstadt 
West Germany 
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TABLE VI -4 

COMPOSITES -RELATED VISITS TO RELEVANT ORGANIZATIONS 
by RPI Faculty/Staff/Students 
(April 30, 1984 through September 30, 1984) 


Visited Date By. Purpose 


East China Insti- 
tute of Textile 
Science and Tech- 
nology, Shanghai, 
China 

5/12- 

23 

Prof. B. 
Wunderlich 

Attended series of ten 
lectures 

Jiao Tong Univer- 
sity, Shanghai, 
China 

5/17 

Prof. B. 
Wunderlich 

Attended seminar 

Fudan University 
Shanghai, China 

5/22 

Prof. B. 
Wunderlich 

Attended lecture 

Aerospace Corpor- 
ation, El Segundo, 
CA 

5/22- 

23 

Prof . R. J . 
Diefendorf 

Attended Rocket Propul- 
sion Laboratory Work- 
shop on "Carbon/Carbon 
Exit Cones" 

Jing-shan Chemical 
Complex, Shanghai, 
China 

5/24 

Prof. B. 
Wunderlich 

Attended seminar 

Goodrich Research 
and Development 
Center, Breckville, 
OH 

5/25 

Prof. R. J. 
Diefendorf 

Attended seminar 

University of Sci- 
ence and Technol- 
ogy of China, 

Hefei , China 

5/26- 

27 

Prof. B. 
Wunderlich 

Attended two lectures 
and seminar 

AVCO Corporation 
Specialty Materi- 
als Division, 
Lowell, MA 

5/30 

Prof. N. S. 
Stolof f 

Attended lecture 

Beijing Institute 
of Technology, 
Beijing, China 

6/2-5 

Prof. B. 
Wunderlich 

Attended three lectures 

General Motors 
Technology Center, 
Warren, MI 

7/25 

Prof. R. J. 
Diefendorf 

Attended seminar 

Lord Corporation, 
Erie, PA 

8/28 

Prof. R. J. 
Diefendorf 

Attended seminar 



TABLE VI-4 continued 


Visited 


General Electric 
Comp any , Water- 
ford, NY 

General Motors 
Technology Center, 
Warren, MI 


Date 

By 

9/13 

Prof. R. J. 


Diefendorf 

9/28 

Prof. R. J. 


Diefendorf 


Purpose 


Attended seminar 


Attended seminar 
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reports, usually on the part of one of the Co-Principal In- 
vestigators. Off-campus visitors, at RPI during a BBL day, 
are often invited to "sit in". Table VI-5 lists a calender 
of internal, oral progress reports as they were given at 
BBL's during this reporting period. 

As indicated in the Introduction of this progress report, 
an initiative has been implemented which is intended to bring 
about increased communication between NASA researchers and 
their RPI counterparts. A first step in that direction has 
been taken by the planning of a series of Research Coordina- 
tion Meetings. The first of such meetings was held at RPI 
on August 27, 1984, and fourteen members of RPI's Composite 
Materials and Structures Program along with seven members of 
NASA Langley Research Center's Materials Division were pres- 
ent. The day-long meeting began with a welcoming statement by 
Dr. R. Loewy, RPI's Composite Program Co-Principal Investi- 
gator, followed by introductory remarks by Dr. L. Vosteen of 
NASA Langley Research Center. The meeting was then split 
into groups of participants with parallel research interests. 
Participants were free to move from group to group as they 
so desired. A final session with all participants provided 
an opportunity for individuals to make comments or suggestions 
with regard to furthering cooperation between organizations. 
The agenda for this meeting is shown in Table VI-7. A second 
of such meetings with researchers from NASA Lewis Research 
Center is in the planning stage. 
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TABLE VI -5 



COMPOSITE MATERIALS AND STRUCTURES PROGRAM 


BROWN BAG LUNCH (BBL) SCHEDULE 



(April 30, 1984 through September 30, 

1984) 

Date 

Topic 

Responsible 

Faculty 

5/4 

Administrative Report 

R. Loewy 


Matrix Characterization and Environ- 
mental Effects 

S. Sternstein 


Fabrication Technology Experiments 

F . Bundy 
H. Hagerup 
V. Paedelt 

9/7 

Administrative Report 

R. Loewy 


Report on NASA Langley Visit 

R. Loewy 

9/14 

Administrative Report 

R. Loewy 


Edge Failures 

S . Sham 


Beams and Warping 

0. Bauchau 

9/28 

Administrative Report 

J. Diefendorf 


Eutectucs 

N. Stoloff 


Fabrication Technology Experiments 

F . Bundy 
H. Hagerup 
V. Paedelt 
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TABLE VI- 6 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
RESEARCH COORDINATION MEETING 
NASA Langley Research Center and Rensselaer Polytechnic Institute 

August 27/ 1984: AGENDA 


1st Plenary Session 
10:45 Welcome 

11:00 Introductory Remarks 

Individual Meetings (flexible) 

11:30-12:30 Groups A, B, C 

12:30-1:30 Lunch 

1:30-4:00 Groups A, B, C, D 

2nd Plenary Session 

4:00 General Discussion 


Dr. R. Loewy 
Rensselaer Polytechnic 
Institute 

Dr. L. Vosteen 

NASA Langley Research 

Center 


5:00 


Adjourn 
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TABLE VI-7 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
RESEARCH COORDINATION MEETING 
NASA Langley Research Center and Rensselaer Polytechnic Institute 

August 27, 1984; ATTENDEES 


Dr. R. Loewy, Rensselaer Polytechnic Institute 
Dr. L. Vosteen, NASA Langley Research Center 

Group A ; Fiber and Matrix Development and Laminate Processing 

NASA - J. Maas 

J. Hinkley 

RPI - J. Diefendorf 

M. Shephard 
S. Sternstein 
B. Wunderlich 

N. Stoloff 

Group B : Fracture and Fatigue 

NASA - W. Elber 

K. O'Brien 

RPI - E. Krempl 

M. Shephard 
S . Sham 

Group C ; Mechanics and Structural Analysis 

NASA - M. Shuart 
M. Nemith 

RPI - D. Goetschel 
M. Shephard 

O. Bauchau 

Group D ; Innovative Fabrication Techniques (Sailplane Projects) 

NASA - M. Shuart 
M. Nemith 

E. Elber 

RPI - D. Goetschel 
V. Paedelt 
0. Bauchau 

F. Bundy 
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PERSONNEL 


Co-Principal Investigators 

Ansell, George S., Ph.D. 
Loewy, Robert G., Ph.D. 
Wiberley, Stephen E., Ph.D. 

Senior Investigators 

Bauchau, 0. , Ph.D. 
(Structural dynamics, ad- 
vanced composites) * 

Bundy, F. P. , Ph.D. 

(Physical chemistry and 
structures testing) * 

Diefendorf 1 , R. J. , Ph.D. 
(Fabrication, resin matrix, 
fiber behavior, interfaces)* 

Feeser 1 , L. J. , Ph.D. 
(Computer applications and 
graphics, computer-aided 
design, optimization)* 

Goetschel, D. B. , Ph.D. 
(Structural analysis design 
and testing) * 

Hagerup, H. J. , Ph.D. 
(Aerodynamics, configura- 
tion, pilot accomodation, 
flight testing) * 

Krempl , E. , Dr.Ing. 

(Fatigue studies, failure 
criteria) * 

Sham, T-L. , Ph.D. 

(Fracture mechanics, com- 
posites) * 


Dean, School of Engineering 
Institute Professor 
Professor of Chemistry 

Assistant Professor of Aero- 
nautical Engineering 


Research Professor of Mate- 
rials Engineering 


Professor of Materials Engi- 
neering 


Professor of Civil Engineering 
Associate Dean, School of Engi- 
neering 


Assistant Professor of Me- 
chanical Engineering 

Associate Professor of Aero- 
neering Engineering 


Professor of Mechanics and 
Director of Cyclic Strain 
Laboratory 

Assistant Professor of Me- 
chanical Engineering 


Fields of Speciality 

^■Member of Budget Committee together with Co-Principal Inves- 
tigators 
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Shephard, M.S., Ph.D. 
(Computer graphics, finite 
element methods)* 

Sternstein 1 , S. S., Ph.D. 
(Failure analysis, matrix 
behavior, moisture effects) * 

Stoloff, N. S. Ph.D. 
(Mechanical behavior of 
crystals, order-disorder re- 
actions, fracture and stress 
(corrosion) * 

Wunderlich, B. , Ph.D. 
(Processing science, consti- 
tuent material characteris- 
tics) * 

Research Staff 


Associate Director, Center for 
Interactive Computer Graphics 
and Assistant Professor of 
Civil Engineering 

William Weightman Walker 
Professor of Polymer Engineer- 
ing 

Professor of Materials Engi- 
neering 


Professor of Chemistry 


Manager & Master Technician, 
Paedelt, Volker 

Research Associates 
Grebowicz, Janusz, Ph.D. 

Research Administrator 
Trainer, Asa, M.S. 

Graduate Assistants 
An Duek , M . S . 

Burd, Gary, M.S. 

Chen, . Kuong- jung, M.S. 
Coffenberry, Brian, B.S. 
Gjaja, Niko, M.S. 

Hu, Tsay-hsin, M.E. 

Judovits, Lawrence, M.S. 

Liu, Shiann-hsing, M.S. 
Shick, Douglas, M.S. 


Composites Laboratory 


Singh, Sachchida, M.S. 
Srinivasan, Krishna, B.Tech. 
Szewczyk, Christine, B.S. 

Uzoh, Cyprian, B.S. 

Weidner, Theodore, B.S. 

Yang, Chung-Ming (Philip) , M.S. 
Yehia, Nabil, M.S. 

Yerry, Mark, M.E. 

Yurgartis, Steven, M.S. 


Fields of Speciality 

1 Meinber of Budget Committee together with Co-Principal Inves- 
tigators 



Undergraduate Assistants 

- Seniors 

Basel, Roger 

Kim, S. Kwong 

Cimino, Paul 

Kirker, Philip 

Di Lello, Frank 

Krupp, Alan 

Father, Richard 

Lowy, Robert 

Hubner, Angela 

Mao , Marlon 

Karkow, Jon 

Williams, Thomas 

Undergraduate Assistants 

- Juniors 

Bell, Joseph 

Kashynski, Stephen 

Burdick , Mark 

O'Connell, James 

Douskay, Eugene 

Ragczewski, David 

Egbert, Mark 

Young, Richard 

Hill, Stephen 


Undergraduate Assistants 

- Sophmores 

Baldwin, Reid 

Park, Brian 

Meyer, John 

Rosario, Estrella 
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